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Abstract 

The  Holloman  High  Speed  Test  Track  (HHSTT)  is  a  rocket-powered  sled  track 
facility  used  for  testing  a  variety  of  hypervelocity  aerospace  applications.  The  current 
speed  record  is  6,453  miles  per  hour.  While  this  seems  fast  there  are  customers  at  the 
track  that  have  requirements  demanding  even  faster  speeds.  Significant  research  has  been 
conducted  in  the  area  of  rail  gouging  as  it  relates  to  the  test  track,  and  efforts  are  under 
way  to  reduce  and  even  eliminate  this  phenomenon.  Any  steps  taken  in  this  effort  may 
eliminate  catastrophic  sled  failures  caused  by  gouging,  however  wear  is  another 
damaging  issue  that  needs  to  be  understood.  This  research  evaluates  wear  in  two 
fashions.  Eirst,  data  from  the  Dynamic  Analysis  and  Design  System  (DADS)  software, 
used  by  the  HHSTT  to  model  sled  loading  and  vibration,  is  evaluated  in  a  theoretical 
model  originating  from  “high-speed”  pin-on-disk  experimentation.  The  second  method 
evaluates  wear  by  conducting  a  series  of  simulations  in  which  temperature,  rail  geometry 
in  the  form  of  an  asperity,  speed  and  rail  coating  are  varied.  These  short  duration 
simulations  are  performed  with  CTH,  a  hydrocode  designed  for  analyzing  hypervelocity 
impact  problems.  A  validated  CTH  model,  used  in  previous  gouging  research,  was  used 
to  estimate  wear  by  examining  local  areas  of  failure  and  melting.  DADS  provides  forces 
related  to  frictional  heating  and  vertical  velocity  information  that  will  be  incorporated  as 
input  into  CTH. 
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I.  Introduction 


Holloman  High  Speed  Test  Track 

The  Holloman  High  Speed  Test  Track  (HHSTT)  located  at  Holloman  Air  Force 
Base  in  Alamogordo,  New  Mexico  is  the  world’s  premier  rocket  sled  test  facility  that 
began  testing  aerospace  applications  in  1950.  Rocket  powered  sleds  are  propelled  down 
a  track  of  almost  51,000  feet  on  heavy-duty  crane  rails  (171  Ib/yrd)  performing  tests 
investigating  hypersonic  environments,  aircraft  ejection  seats  and  munitions  and 
aerodynamic  related  effects.  The  track  provides  cost,  safety,  efficiency  and  observability 
advantages  over  flight-testing  while  it  also  offers  a  more  realistic  approach  to  laboratory- 
based  tests.  The  fastest  sled  on  record  reached  a  speed  of  6,453  miles  per  hour  in  April  of 
2003  on  a  newly  developed  dual  rail  sled  system.  As  technology  and  interest  expand, 
customers  are  expressing  interest  in  achieving  even  faster  speeds.  In  an  effort  to  provide 
consistent,  reliable  and  payload  safe  environments,  the  track  engineers  are  forced  to  study 
and  understand  the  effects  of  these  hypersonic  speeds. 

The  sled  is  connected  to  the  rails  by  four  “shoes”  or  “slippers”  that  wrap  around 
the  railhead  and  act  as  a  guide  for  the  sled  (Figure  1)  as  it  is  propelled  down  the  track. 

The  rails  are  constructed  of  1080  steel.  A  typical  shoe  is  made  of  a  high  strength  steel 
alloy,  VascoMax  300,  and  is  dimensionally  20.32  cm  (8  inches)  in  length,  by  10.16  cm  (4 


1 


inches)  wide  and  2.54  cm  (1  in)  in  thickness.  The  shoes  are  designed  to  have  a  nominal 
0.3175  cm  (0.125  in)  gap  between  the  shoe  and  the  rail.  Sleds  are  eonstrueted  with  one 
or  more  roeket  motors  depending  on  the  size  of  the  payload  and  speed  requirements. 
Frequently  pusher  sleds,  earrying  additional  motors  that  aet  in  stages  to  propel  the  sled 
system  down  the  traek,  are  used. 


As  the  sled  veloeity  inereases,  aerodynamic  forces  and  heating  beeome 
significant.  At  supersonic  speeds,  the  leading  edges  of  the  shoes  are  exposed  to  elevated 
temperatures  as  they  are  behind  a  strong  shook  wave.  Aerodynamio  lift  assists  in  oausing 
the  sled  to  osoillate  or  “bounoe”  as  it  moves  down  the  traek.  In  effeot,  the  sled  is 
oonstantly  transitioning  from  making  oontaot  with  the  top  surfaoe  of  the  rail  to  the  shoe 
keeping  it  from  flying  off  of  the  rail  by  oontacting  the  underside  of  the  rail.  This 
“bounoing”  is  not  only  oaused  by  aerodynamio  foroes  but  also  by  the  rail  irregularities. 
These  irregularities  are  a  result  of  inherent  manufaoturing  and  alignment  toleranoes  of  the 
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rails.  As  a  result  the  gap  between  the  shoe  and  the  rail  is  eonstantly  changing.  As  the 
sled  speed  increases,  this  changing  gap  causes  severe  vibrations  to  occur  that  only 
intensify  with  speed  [12:6].  In  an  effort  to  reduce  the  harshness  of  the  environment  so 
that  higher  speeds  could  be  attainable,  a  new  dual  rail  system  was  developed  under  the 
Hypersonic  Upgrade  Program  (HUP)  [11].  Figure  2  shows  the  configuration  of  a  dual 
rail  sled  of  which  this  research  focuses. 


Figure  2  -  Dual  Rail  Sled  with  Pusher  Motors 


Background 

The  Air  Force  Institute  of  Technology  (AFIT)  has  been  studying  the  HHSTT  for 
several  years,  focusing  all  of  the  efforts  on  rail  gouging.  Gouging  can  be  characterized  as 
mixing  in  a  fluid-like  behavior  of  two  materials  under  sliding  contact  that  experience 
plastic  deformation  at  high  strain  rates.  Laird  provides  a  complete  definition  of  gouging 
in  his  work  [17:6].  Cinnamon’s  research  on  gouging  concluded  that  “a  new  rail 
alignment  criteria... should  mitigate  future  incidences  of  hypervelocity  gouging”  [7:10- 
11].  If  this  proves  to  be  true,  and  gouging  is  eliminated,  the  next  problem  that  is  of 
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concern  to  the  engineers  at  Holloman  is  that  of  wear.  In  the  problem  of  wear,  the  focus 
changes  from  looking  at  the  rail  to  examining  the  behavior  of  the  shoe  as  it  slides  along 
the  rail. 

Wear 


In  1976,  it  was  estimated  that  over  $24  billion  was  spent  on  replacing  worn  items 
on  automobiles  alone  [21:7-8].  Many  things  in  addition  to  automobiles  are  affeeted  by 
the  wear  phenomenon.  Wear  is  not  only  a  financial  matter;  it  often  raises  safety  concerns 
in  sueh  applications  as  vehicular  brakes.  Although  wear  is  generally  seen  in  a  negative 
light,  it  also  has  some  useful  purposes  sueh  as  grinding,  sanding  and  polishing  or  even  the 
simple  exercise  of  writing  with  a  pencil.  The  phenomenon  of  wear  warrants  a 
comprehensive  understanding;  unfortunately  wear  is  very  eomplex  and  multifaceted. 

Wear  fits  under  a  larger  umbrella  of  science  ealled  tribology.  Generally  speaking, 
tribology  is  the  study  of  the  behavior  of  materials  and  associated  particles  in  contact  as 
they  undergo  relative  motion.  Friction  and  lubrication,  along  with  wear,  are  the  major 
fields  included  in  study  of  tribology.  Wear,  aecording  to  I.M.  Hutchings,  is  defined  as 
“the  loss  of  material  from  a  surfaee,  transfer  of  material  from  one  surfaee  to  another  or 
movement  of  material  within  a  single  surfaee”  [13:1].  For  the  purposes  of  this  researeh, 
wear  will  be  defined  as  the  loss  of  material  due  to  failure  and  or  melting. 

Signifieant  time,  money  and  effort  have  been  spent  on  the  study  of  wear,  and  to 
this  day  it  is  still  not  completely  understood.  Bowden  and  Tabor  summarize  the 
mystifying  behavior  of  wear  in  the  following  paragraph: 
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“There  are  no  simple  laws  of  wear  as  there  are  for  frietion.  On  the  whole, 
hard  solids  wear  less  than  soft,  though  polyethylene  wears  less  than  steel. 

The  wear  generally  inereases  with  the  distanee  rubbed,  though  this  may 
not  be  true  if  appreeiable  baek  transfer  oeeurs.  The  wear  generally 
inereases  with  load,  though  this  may  not  be  true  if  more  severe  running 
eonditions  produee  appropriate  struetural  ehanges  in  the  surfaee  layer.”  - 
[5:148] 

The  eomplexity  of  wear  may  be  attributed  to  the  ever-ehanging  environment  in 
whieh  two  materials  interaet.  The  proeess  of  wear  ehanges  the  surfaee  geometry  of  the 
materials.  Inereasing  temperatures  due  to  frietional  heating  and  plastieity  eause  material 
strength  properties  and  hardness  to  ehange,  and  if  enough  heat  is  generated,  a  material 
phase  ehange  may  oeeur.  This  dynamie  environment  signifieantly  impaets  the  type  and 
amount  of  wear  that  oeeurs.  Beeause  of  this,  no  two  systems  are  identieal.  As  Kato 
suggests,  “wear  is  not  a  material  property,  it  is  a  systems  response.”  [16:9] 

Wear  ean  be  eharaeterized  as  a  meehanieal,  ehemieal  or  thermal  proeess. 
Meehanieal  wear,  defined  as  deformation  and/or  fraeturing  of  the  surfaee,  ean  be 
represented  by  abrasion,  adhesion,  fatigue  and  plastie  flow.  Chemieal  wear,  also  referred 
to  as  triboehemieal  or  eorrosive  wear,  is  eharaeterized  by  a  growing  “ehemieal  reaetion 
film”  whieh  is  “aeeelerated  meehanieally  by  frietion”  [16:1 1].  Thermal  wear  is  best 
represented  by  surfaee  melting  eaused  by  heat  generated  from  frietion,  however  brittle 
material  wear  resulting  from  thermal  shoeking  is  also  ineluded  in  thermal  wear  [16:1 1]. 

The  following  paragraphs  inelude  a  brief  diseussion  on  the  eoneepts  of  wear  on 
the  mierostruetural  level  and  how  it  is  affeeted  by  the  environment  in  whieh  it  takes 
plaee.  This  diseussion  is  entirely  relevant  to  sliding  speeds  below  100  m/s.  It  is  beeause 
wear  is  less  eommon  and  eertainly  less  understood  at  speeds  higher  than  100  m/s  that 
warrants  this  researeh.  Therefore  the  eonelusion  of  this  ehapter  provides  an  engineering 


5 


approach  for  determining  wear  for  speeds  that  are  of  interest  to  the  HHSTT,  that  is  speeds 
on  the  order  of  3000  m/s. 

Adhesion  -  Dry  Sliding  Wear 

Dry  sliding  wear  consists  of  two  surfaces  moving  relative  to  each  other  without 
the  presence  of  particles,  either  loose  or  attached  to  one  or  both  of  the  surfaces  and 
without  a  lubricant  film.  This  type  of  wear  can  be  thought  of  as  adhesion  because  it  is 
heavily  dependent  on  the  attraction  of  the  mating  materials.  Abrasive  wear  on  the  other 
hand  is  classified  as  sliding  wear  with  particles,  which  will  be  discussed  briefly  in  the 
next  section. 

As  two  materials  come  together,  only  the  peaks  of  their  surface  asperities  contact, 
resulting  in  highly  loealized  stresses.  Figure  3  is  a  depiction  of  a  ground  steel  surfaee, 
shown  here  to  demonstrate  the  level  of  contact  that  two  similar  mating  surfaces  undergo. 


Figure  3  -  Typical  Surface  Illustrating  Asperities  [4:149] 

Wear  comes  about  as  pieces  of  these  asperities  are  deformed  or  removed.  The 
mechanism  that  brings  about  the  eventual  particle  removal  is  not  completely  understood. 
Many  theories  have  been  suggested;  all  of  them  are  foeused  on  the  combination  of 
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normal  stress  bought  about  by  loading;  shear  stress  brought  about  by  transverse  motion 
and  friction;  and  temperature  that  increases  with  both  load  and  velocity.  Some  theories 
suggest  that  particles  break  off  as  a  result  of  plastic,  plate-like  movement  ending  in 
fracture  of  the  attached  end  as  in  Figure  4a,  while  others  suggest  that  layers  fracture 
independently  due  to  shear  instead  of  breaking  off  in  a  stack  (Figure  4b). 


Figure  4  -  Adhesive  Wear  -  Plasticity  Resulting  in  Plate-Like  Debris  [14:99] 

Some  propose  that  smaller  wear  particles  transfer  from  one  surface  to  another  from 
adhesion.  These  particles  build  up  over  time  and  eventually  separate  as  debris  as  in 
Figure  5. 


(a)  Initial  asperity  contact  (b)  Junction  failure 


(c)  Transfer  of  particle  (d)  Secondary  contact 


(e)  Merging  of  two  particles  (f)  Transfer  particle  just 

before  detachment 


Figure  5  -  Smaller  Particle  Debris  Form  as  a  Result  of  Adhesion  [14:100] 
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There  are  also  those  that  propose  eracks  generate  below  and  parallel  to  the  surface 
eventually  resulting  in  the  delamination  of  a  particle.  All  of  these  theories  deal  with  the 
asperities  reaching  the  yield  point,  deforming  plastically  and  eventually  fracturing. 
Examination  of  worn  materials  has  confirmed  that  transfer  layers,  laminate  debris  and 
subsurface  cracking  exist,  validating  all  of  the  theories,  however  the  exact  process 
remains  uncertain  [14:92-103]. 


Abrasion 

The  previous  section  focused  on  wear  without  particles.  Abrasion  occurs  in  one 
of  two  forms:  particles  attached  on  one  of  the  mating  surfaces  (Figure  6a)  or  particles 
existing  unattached  between  the  two  mating  surfaces  (Figure  6b). 


(b) 


Figure  6  -  Two  Forms  of  Abrasion  [14:133] 

If  wear  debris  is  generated  from  adhesion,  as  described  above,  and  not  ejected 
from  the  interface,  the  result  is  abrasion.  The  existence  of  abrasive  particles  in  the  form 
of  wear  debris  has  been  shown  to  increase  the  coefficient  of  friction  between  the  two 
surfaces.  The  coefficient  of  friction  is  a  function  of  the  debris  particle  size,  which  is 
dependent  on  the  hardness  of  the  materials  [9:137].  The  mode  of  abrasion  is  mostly 
dependent  on  the  shape  of  the  abrasive  particle. 
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Figure  7  -  Three  Modes  of  Abrasion  [14:147] 

Figure  7  depicts  the  three  modes  as  cutting  (a),  ploughing  (c)  ,  and  the  intermediate  mode 
called  wedge  formation  (b).  Cutting  is  analogous  to  that  of  a  single  point  cutting  tool 
which  results  in  the  removal  of  a  chip.  Ploughing  is  like  a  rolling  wave,  however  this 
mode  does  not  result  in  debris  particles  whereas  the  wedge-forming  mode  does.  These 
modes  mostly  depend  on  the  face  angle,  0,  of  the  abrasive  particle  [14:147-148]. 
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Melt 


This  mechanism  of  wear  is  encountered  when  the  interface  temperature  reaches 
that  equal  to  one  or  both  of  the  melting  temperatures  of  the  mating  materials.  This  is  a 
result  of  the  frictional  heat  generated  from  sliding  which  is  a  function  of  load  and 
velocity.  It  is  generally  accepted,  in  theory,  that  as  a  molten  layer  develops,  the 
coefficient  of  friction  decreases  due  to  the  hydrodynamic  layer,  similar  to  lubrication. 

For  metals  with  low  melting  points  like  bismuth,  experimental  evidence  exists  for  the 
wear  rate  to  increase  when  melting  occurs  [23:610],  however  melt  wear  models  are  not 
conclusive  across  the  board. 

Oxidation 

This  mechanism  of  wear  has  to  do  with  the  rapid  build-up  of  oxide  layers.  During 
the  process  of  oxidation,  a  chemical  reaction  occurs  between  the  base  metal  and  the 
surrounding  oxygen,  resulting  in  the  formation  of  oxides  of  that  base  metal.  Oxides  in 
general,  are  harder,  more  brittle  and  have  a  lower  shear  strength  than  the  base  metal. 
These  considerations  can  increase  the  wear  rate  if  multiple  passes  occur.  These  oxide 
layers  are  built  up  and  then  broken  down  when  they  reach  a  critical  thickness.  Lower 
temperatures  from  slower  sliding  speeds  cause  oxides  to  build  up  only  on  the  asperities, 
whereas  higher  temperatures  from  higher  speeds  and  loads  cause  a  more  general  oxide 
layer  to  form  [14:103-105]. 
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Wear  Maps 


As  materials  experience  repeated  contact  during  the  wearing  process,  the 
microstructure,  micro-surface  geometry  and  material  chemistry  change  as  seen  in  the 
previous  sections.  These  changes  bring  about  different  wear  mechanisms  that  dominate 
during  the  wearing  process.  In  order  to  track  this  dynamic  process  in  terms  of  wear,  the 
idea  of  a  wear  a  map  was  developed.  A  wear  map  is  a  graphical  representation  used  to 
indicate  operational  boundaries  such  as  load  and  speed  for  different  wear  mechanisms 
and/or  rates.  Bayer  [2:383]  also  mentions  a  transition  diagram  as  a  form  of  wear  map. 
These  diagrams  are  generated  from  empirical  data  by  plotting  experimental  wear 
quantities  against  two  operating  parameters  to  designate  different  wear  regimes.  A  wear 
map  or  transition  diagram  for  steel  in  a  pin-on-disk  configuration  can  be  seen  in  Figure  8 
as  a  function  of  normalized  pressure  and  normalized  velocity.  These  quantities  will  be 
defined  in  the  next  chapter  (see  equations  (5)  and  (6)). 
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Figure  8  -  Wear  Map  for  Metals  [18:21] 


In  the  attempt  to  determine  the  total  wear  of  a  surface,  it  will  be  necessary  to  account  for 


the  wear  from  all  of  the  mechanisms. 


Engineering  Solution  to  Wear 

As  previously  discussed,  experimental  wear  maps  show  where  certain  wear 
mechanisms  will  dominate  in  terms  of  some  operational  parameters.  It  is  therefore 
desirable,  in  the  absence  of  experimental  data,  to  be  able  predict  and  quantify  wear  and 
its  mechanisms  for  a  given  wear  scenario.  This  can  be  done  by  analytically  developing  a 
wear  model  from  a  specific  system  and  applying  it  generically  to  others. 

Wear  is  generally  expressed  in  the  form  of  a  rate,  volume  of  material  removed  per 
distance  slid.  The  wear  rate,  as  alluded  to  in  the  last  section,  does  not  remain  constant 
throughout  the  life  of  the  system.  Typically  there  is  a  “break-in”  period  where  the  wear 
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rate  is  high,  and  then  it  tapers  off  with  time  as  the  mating  surfaces  are  smoothed.  In  the 
case  of  the  Holloman  sled,  the  wear  rate  may  never  reach  a  steady  state  value  due  to  the 
constant  acceleration  and  intermittent  contact  with  the  rail.  The  volume  of  wear  can  be 
expressed  as 


V  =  k,'^ 
H 


(1) 


where  H  is  the  material  hardness  in  units  of  pressure,  Fff  is  the  normal  load,  d  is  the 
sliding  distance  and  kA  is  the  dimensionless  wear  coefficient.  The  wear  rate  can  be 
expressed  as 


V 

w  =  — 

d 


(2) 


or  as  a  dimensionless  quantity,  w 


w  = 


V 

And 


(3) 


with  An  being  the  nominal  area  of  contact.  Substituting  equation  (1)  into  equation  (2) 
results  in  the  most  rudimentary  form  of  a  wear  model. 

w  =  k^^  (4) 

H 

This  equation  is  known  as  Archard’s  wear  equation,  and  kA  is  known  as  Archard’s  wear 
coefficient,  named  after  a  British  scientist  who  was  one  of  the  first  to  derive  the 
relationship  between  wear  and  load  in  1953.  Archard  developed  this  equation  based  on 
the  probability  of  an  asperity  in  contact  producing  a  wear  particle.  He  took  no  specific 
consideration  for  the  wear  mechanism  [14:83-84].  Others  have  derived  a  wear  equation 
identical  to  equation  (4)  by  beginning  with  the  assumptions  of  adhesive  wear  as  well  as 
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abrasive  wear  [16:1 1-12],  Considering  that  almost  all  wear  seenarios  involve  these  two 
meehanisms,  equation  (4)  ean  be  generally  applied  as  an  estimation  of  wear.  In  that  light, 
kA  is  left  to  take  on  a  large  range  of  values  depending  on  the  type  of  wear  under 
eonsideration. 

Problem  Statement  and  Objectives 

As  the  Holloman  High  Speed  Test  Traek  eontinues  to  push  the  speed  boundaries, 
the  problem  of  wear  beeomes  a  significant  issue.  Wear  in  general  is  complex  and  not 
completely  understood;  even  less  is  known  about  wear  at  speeds  above  500  m/s.  With 
Holloman  hoping  to  soon  achieve  3,000  m/s,  this  problem  becomes  uniquely  interesting. 
The  HHSTT  hopes  to  be  able  to  predict  the  shoe  wear  so  that  it  does  not  completely  wear 
through  during  a  test  and  gouge  the  rail  or  worse,  cause  a  catastrophic  sled  failure,  as  the 
shoe  would  loose  its  ability  to  keep  the  sled  on  the  track.  Sled  failures  are  not  only  costly 
to  the  track  in  terms  of  maintenance  and  repair;  they  are  costly  to  the  customer  as  well. 
The  objectives  of  this  research  are  to: 

•  Relate  low  speed  wear  experiments  to  higher  velocity  applications  using 
parameters  such  as  pressure  and  velocity. 

•  Evaluate  high  velocity  wear  using  CTH  considering  the  effects  of  temperature, 
rail  geometry,  velocity  and  coatings. 

Methodology  Preview 

This  research  consists  of  the  development  of  two  models,  an  analytical  model  and 
a  numerical  model.  The  numerical  model  is  created  in  CTH  while  the  analytical  model  is 
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developed  from  the  work  of  Lim  and  Ashby  [18]  and  applied  to  the  data  provided  by  the 
Dynamie  Analysis  and  Design  System  (DADS). 

CTH,  a  hydroeode  that  was  developed  by  Sandia  National  Labs  for  modeling 
high-energy  impact,  has  been  used  to  successfully  model  the  problem  of  rail  gouging  at 
Holloman.  CTH  was  originally  chosen  to  model  gouging  because  of  its  ability  to  model 
high  pressure  interaction  between  materials  leading  to  large  deformations  and  the 
propagation  of  shock  waves.  This  was  also  used  to  estimate  shoe  wear  considering  the 
sled  contacting  the  rail  and  generating  areas  of  higher  pressure.  While  gouging  is 
characterized  by  the  flow  and  mixing  of  the  shoe  and  rail  materials,  wear  is  characterized 
by  localized  melting,  failure  and  debris  generation.  CTH  will  be  used  to  quantify  wear 
by  looking  at  the  temperature  and  failure  zones  in  the  shoe,  as  well  as  material  loss,  if  any 
exists. 

Computer  generated  sled  parameters  such  as  force  and  velocity  as  functions  of 
time  are  available  through  a  sled  as  modeled  in  DADS.  The  particular  DADS  file  that 
was  used  for  this  research  was  a  dual  rail  sled  with  a  mass  of  841  kg  (1854  lb);  it  is 
partially  listed  in  Appendix  D.  The  model  was  based  on  an  actual  sled  run  from  October 
2000,  designated  Mission  80X-A1,  but  was  modified  to  travel  farther  and  faster.  The 
model  traveled  a  distance  of  over  3.65  km  (12,000  ft)  reaching  a  maximum  speed  of 
approximately  3  km/s  (10,000  ft/s)  in  2.5  seconds.  A  wear  rate  can  be  generated  for  each 
contact  the  sled  makes  with  the  rail  and  then  integrated  over  the  distance  to  determine  the 
wear  volume  and  therefore  thickness.  The  DADS  also  provides  a  method  for  calculating 
frictional  heating  that  will  be  used  to  calculate  a  temperature  distribution  in  the  shoe. 
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CTH  is  used  to  run  a  series  of  simulations  at  eertain  veloeities  to  determine  the 


effeets  that  temperature,  rail  eoatings,  and  rail  diseontinuities  in  the  form  of  an  asperities, 
have  on  wear.  The  temperature  input  is  the  result  of  the  DADS  heat  transfer  solution  just 
mentioned.  Beeause  CTH  runs  for  very  short  time  intervals,  on  the  order  of  10 
mieroseeonds,  only  seleeted  seenarios  ean  be  examined.  The  goal  is  to  eventually 
eharaeterize  the  entire  sled  run  by  a  broad  range  of  seenarios  and  statistieally  ealeulate 
the  total  wear.  The  two  methods  for  evaluating  wear  ean  be  used  to  eheek  eaeh  other, 
however  it  will  not  be  until  aetual  test  traek  wear  data  is  gathered  for  eomparison  that  the 
two  methods  will  be  validated. 
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II,  Methodology 


Pin-On-Disk 

Pin-on-disk  experiments  have  been  used  for  many  years  to  simulate  wear  in  the 
laboratory.  These  tests  can  vary  in  geometric  setup;  for  example,  the  pin  can  be  applied 
to  the  edge  of  the  disk  instead  of  the  face.  Bayer  [2]  dedicates  several  chapters  to  the 
discussion  of  wear  testing  and  techniques.  This  research  focuses  on  the  setup  illustrated 
in  Figure  9.  A  pin  of  known  material  properties  and  geometry,  namely  density,  hardness, 
length,  and  diameter,  is  placed  against  a  disk  with  constant  angular  velocity  for  some 
defined  period.  From  that  period  and  angular  velocity,  the  sliding  distance  can  be 
determined.  When  simulating  a  linear  sliding  scenario  in  one  direction,  such  as  the 
Holloman  sled,  or  a  projectile  traveling  down  a  gun  barrel,  the  pin  can  be  moved  in  a 
radial  direction  away  from  the  center  of  the  disk  to  allow  it  not  to  pass  over  the  same  path 
on  the  disk  more  than  once.  This  radial  rate  would  also  be  used  to  calculate  the  sliding 
distance.  The  force  applied  to  the  pin  in  the  normal  direction  is  measured  as  well  as  the 
frictional  force  exerted  on  the  pin.  The  frictional  force  can  be  determined  by  calculating 
the  bending  moment  in  the  pin  from  a  strain  gage  located  on  the  front  or  back  of  the  pin. 


Figure  9  -  Pin-on-Disk  Experimental  Setup  [14:79] 
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Sled  vs,  Pin-on-Disk 


The  speed  of  the  sled  alone  makes  the  HHSTT  unique.  There  are  no  other  tests  or 
experiments  that  ean  be  used  as  an  equivalent  comparison  when  considering  only 
velocity.  The  fastest  experimental  data  available  reaches  only  549  m/s  (1800  ft/s).  On 
the  other  hand,  the  pin-on-disk  experimentation  can  achieve  incredible  pressures,  78  MPa 
(1 1,000  psi)  for  example,  by  only  applying  a  246  N  (55  Ibf)  load  to  a  2  mm  (0.08  in) 
diameter  pin.  If  one  would  consider  the  shoe,  having  an  area  of  206.5  cm  (32  in  ),  large 
pressures  are  not  achievable  because  of  the  large  contact  area.  It  would  require  a  load  of 
1.57  X  10^  N  (352  kips)  to  produce  the  same  pressure  as  the  pin-on-disk.  If  one  considers 
the  product  of  pressure  and  velocity,  an  adequate  comparison  can  be  made  between  the 
pin-on-disk  and  the  shoe/rail  system.  A  disk  spinning  at  a  constant  rate  of  45.7  m/s  (150 
ft/s)  with  the  same  246  N  force  applied  to  the  pin  will  produce  a  Pv  value  of  3.57  GPa 
m/s  (1.7  X  10^  psi  ft/s).  In  order  to  achieve  the  same  Pv  value  for  the  sled/shoe  system,  a 
load  of  48.3  kN  (10,860  Ibf)  would  need  to  be  applied  on  the  shoe  traveling  at  a  velocity 
of  1,524  m/s  (5000  ft/s).  For  this  case,  both  the  load  and  the  speed  are  easily  achieved  by 
the  HHSTT.  The  Pv  quantity  appears  frequently  when  studying  wear  because  normal 
force,  surface  area  and  velocity  contribute  to  the  wear  behavior.  As  will  be  shown  later, 
the  range  of  Pv  for  the  pin-on-disk  is  0  to  70  GPa  m/s  and  that  of  the  shoe/rail  is  0  to  35 
GPa  m/s.  Experimental  pin-on-disk  testing  and  models  have  been  sufficiently  developed 
such  that  their  application  to  the  Holloman  shoe,  by  way  of  DADS,  may  be  valid 
considering  initial  assumptions.  A  list  of  these  assumptions  follows: 

•  For  simplification,  one  corner  of  the  four-shoe  sled  will  be  considered.  The  mass 
applied  on  each  corner  is  assumed  to  be  one  quarter  of  the  total  mass  of  the  sled. 
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throughout  the  entire  run.  This  assumption  eliminates  the  effeets  of  pitehing  and 
rolling  of  the  sled. 

•  The  eontaet  between  the  shoe  and  rail  is  uniformly  distributed  over  the  entire  shoe 
area,  206.5  em^  (32  in^). 

•  Contaet  other  than  the  top  of  the  rail  with  the  bottom  of  the  shoe  is  irrelevant  to 
this  researeh. 

•  The  understanding  of  wear  at  low  speeds  ean  be  applied  to  wear  at  high  speeds. 

Lim  and  Ashby  [18] 

In  1987  two  researehers,  Lim  and  Ashby,  eontending  that  the  meehanisms  of  wear 
ehange  during  the  wearing  proeess,  attempted  to  define  a  wear  map  in  order  to  traek  and 
eharaeterize  the  ehanging  wear  meehanisms  and  ultimately  eompute  the  total  wear.  An 
example  of  a  Lim  and  Ashby  wear  map  was  shown  in  Figure  8,  where  the  different  wear 
meehanisms  are  displayed  in  regions  as  funetions  of  normalized  pressure  and  veloeity. 
Their  work  was  based  on  experimental  data  that  was  eolleeted  from  various  sources. 

Most  of  the  data  resulted  from  pin-on-disk  experimentation  such  as  that  published  by 
Montgomery  in  1975  [19].  With  that  in  mind,  it  is  important  to  note  that  Lim  and  Ashby 
developed  their  wear  models  around  the  pin-on-disk  geometry  of  Figure  9. 

Lim  and  Ashby  almost  exclusively  expressed  wear  in  terms  of  normalized 
quantities  as  defined  earlier  in  equation  (3).  It  is  important  to  note  that  equation  (3)  is 
non-dimensional  but  has  the  units  of  volume  of  wear  per  area  of  contact,  per  distance 
slid,  which  becomes  important  when  trying  to  determine  a  wear  thickness  or  volume. 
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They  also  used  normalized  foree  and  normalized  veloeity  terms.  Normalized  foree  is 
defined  as 


and  normalized  veloeity  as 


F  = 


AH 


(5) 


a 

where  v  is  the  velocity,  To  is  the  radius  of  the  contact  area  (for  the  pin-on-disk  scenario) 
and  a  (m^/s)  is  the  thermal  diffusivity. 

By  plotting  the  experimental  data,  they  were  able  to  surmise  that  the  data 
naturally  fell  into  categories  of  wear  mechanisms.  They  then  formulated  a  theoretical 
model  for  each  mechanism  to  match  the  experimental  data.  They  suggest  that  wear 
occurs  in  three  main  mechanisms,  melt,  oxidation-dominated,  and  plasticity-dominated 
wear  [18:3].  This  research  focuses  on  their  melt  and  plastic  wear  models. 


Melt  Wear 

The  model  for  melt  wear  was  based  on  part  of  the  energy  generated  by  friction 
being  used  to  heat  the  entire  pin  while  another  part  is  focused  on  melting  the  pin  at  the 
interface.  This  fraction  is  defined  by  the  variable  a.  Lim  and  Ashby  used  an  expression 
for  a  based  on  transient  conditions  in  the  disk.  Assumptions  are  made  later  for  a  to  be 
constant  at  14.  The  flux,  q”,  in  W/m^  due  to  frictional  heating  is  defined  as 


q  = 


A. 


(V) 
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where  ju  is  the  coeffieient  of  friction.  Some  of  that  flux  goes  into  heating  the  pin  and 
some  to  melting  the  pin,  as 


dT  V 

aq"=-K  —  +  L^  (8) 

^  dy 

The  first  term  is  taken  from  the  heat  conduction  equation,  which  will  be  examined  later. 

K  is  the  thermal  conductivity  (W/m  K)  and  dT/dy  is  the  temperature  gradient  into  the 
shoe.  The  second  term  is  for  the  portion  of  the  heat  that  melts  the  pin,  where  L  is  the 
latent  heat  of  melting  (J/m  )  and  Vm  is  the  volume  rate  of  melting  (m  /s).  Substituting 
equation  (7)  into  equation  (8)  gives 

afjF^v  _  KjT^-To)  ^  ^  F, 

h  A 

where  4  is  the  equivalent  linear  diffusion  distance.  Lim  and  Ashby  assume  that  all  of  the 
melted  metal  is  removed  as  wear,  therefore  the  equation  for  melt  wear  rate  is 


w  = 


V 

m 


(10) 


Equation  (9)  can  be  solved  for  Vm  and  substituting  it  into  equation  (10)  results  in  the 
following  equation  for  melt  wear  [18:12]. 


w  = 


f  T'  T  ^ 

H  1 

1  A  ) 

ajuFv  ^  ^  -1 

{Tm-To) 


(11) 


where 


Tm  =  melt  temperature 
To  =  sink  temperature 

T*  =  aH  /  K ;  equivalent  temperature  of  the  metal 
P  =  lb/ ro;  dimensionless  parameter  for  bulk  heating 
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Plastic  Wear  /  Mechanical  Wear 


As  mentioned  earlier,  the  most  widely  known  wear  model  is  Arehard’s  law  in  the 
form  of  equation  (4).  By  eombining  equations  (1),  (3)  and  (5),  Arehard’s  law  ean  be 
written  in  dimensionless  form  as 

w  =  k,F  (12) 

Lim  and  Ashby  arrived  at  this  solution  for  what  they  term  plastie  wear  by  suggesting  that 
voids  nueleating  around  material  inelusions  beneath  the  surfaee  generate  wear  partieles  as 
plastie  shear  strain  aeeumulates  beyond  a  eritieal  value  y*  for  speeds  below  0.1  m/s.  In 
this  ease  Ica  is  defined  by 

(13) 

J  A 

with  Yo  being  the  plastie  shear  strain  generated  eaeh  time  that  an  asperity  undergoes 
deformation,/,  is  the  volume  fraction  of  inclusions  and/*  is  the  critical  area  fraction  of 
voids.  By  taking  yo  to  be  1%,/  to  be  10'^  and/*  to  be  0.5,  kA  =  10'^,  a  value 

consistent  with  experimental  results.  Lim  and  Ashby  distinguish  between  mild  and 
severe  plastic  wear  by  changing  /  from  5x10'  to  5  x  10'  respectively,  when  the 
normalized  load,  F  ,  equals  3  xlO'"^  [18:18-19].  Authors  such  as  Hutchings  and 
Rabinowicz  however,  do  not  explicitly  limit  equation  (12)  to  plastic  wear.  They  suggest 
that  this  equation  can  be  applied  to  wear  in  general  [21:155-167;  13:4-5].  For  this  reason, 
wear  calculated  from  equation  (12)  in  this  research  will  be  referred  to  as  mechanical 
wear. 
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Oxidation  Dominated  Wear 


The  equation  for  mild  oxidation-dominated  wear  as  presented  by  Lim  and  Ashby 
is  shown  below  in  equation  (14).  The  development  of  this  equation  is  beyond  the  scope 
but  is  discussed  in  Lim  and  Ashby’s  work  [18:14], 

00 


w  = 


V  y 


exp 


RT, 


F 

V 


(14) 


where 

C  =  volume  oxidized  iron  per  mass  gained  (m  /kg) 

Ao  =  Arrhenius  constant  for  oxidation  (kg  /m  s) 

Zc  =  critical  thickness  for  oxide  film 

Qo  =  activation  energy  for  oxidation  (J/mol) 

R  =  gas  constant  (8.314  J/mol  K) 

Tf  =  flash  temperature 

With  some  basic  information  about  a  wear  scenario,  namely  velocity,  pressure, 
coefficient  of  friction,  geometric  configuration  and  material  properties,  these  models  can 
be  aptly  applied.  However,  due  to  their  complexity,  they  were  first  applied  to  the  data 
published  by  Montgomery  in  1975,  which  listed  experimentally  determined  wear  rates. 
This  way,  the  calculated  values  in  the  above  equations  should  be  close  to  the 
experimental  data.  This  allowed  for  validation  before  using  the  DADS  sled  data. 


Montgomery  Report  [19] 

In  1975  Montgomery  published  a  paper  on  simulated  wear  inside  a  cannon 
muzzle  by  using  the  pin-on-disk  experiment.  The  tests  were  conducted  by  the  U.S.  Army 
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in  the  1940s  and  50s  but  were  elassified  at  the  time.  Montgomery’s  paper  emphasized 
the  use  of  the  produet  of  pressure  and  veloeity,  or  Pv,  for  use  in  analyzing  the  wear 
results.  The  highest  speed  tested  in  the  Army’s  work  was  only  549  m/s  (1,800  ft/s), 
signifieantly  lower  than  what  is  seen  at  Holloman.  However  the  signifieanee  of  the 
Montgomery  data  is  that  it  provides  53  unique  reeords  of  data  eharaeterized  by  veloeity, 
pressure  and  sliding  distanee  with  eorresponding  wear  values.  When  eonsidering  the 
Holloman  sled,  the  eontacts  of  the  sled  with  the  rail  never  oeeur  at  the  same  speed  and 
never  slide  for  the  same  distanee.  This  faet  makes  Montgomery’s  publieation  very 
relevant  due  to  its  wide  variety  of  speed,  load  and  distanee  slid.  Based  on  the  previous 
discussion  of  the  use  of  the  Pv  quantity,  Montgomery’s  data  will  be  used  in  this  fashion 
to  qualify  the  models  that  will  be  applied  to  the  DADS  data. 

Montgomery  concluded  that  the  wear  seen  in  the  Army’s  study  was  categorized  as 
melt  wear  because  the  materials  tested  had  different  melt  temperatures,  and  the  materials 
with  the  lower  melt  temperatures  showed  greater  wear  amounts  [19:297-298].  Lim  and 
Ashby  therefore  used  the  Montgomery  data  to  validate  their  melt  wear  model.  This 
model  will  be  shown  later  in  the  results  section  in  Figure  27. 

DADS 

The  Dynamic  Analysis  and  Design  System,  or  DADS,  is  a  commercial-off-the- 
shelf  software  package  used  to  solve  the  equations  of  motion  for  multi-body  mechanisms. 
The  HHSTT  purchased  DADS  to  help  predict  sled  loading.  Understanding  the  vibrations 
experienced  by  a  sled  was  necessary  to  reduce  sled  failures,  and  design  optimal  sleds  that 
could  protect  sensitive  payloads  at  higher  speeds.  DADS  had  to  be  modified  in  order  to 
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model  the  HHSTT  environment.  In  addition  to  the  aerodynamie  effeets  seen  in  Figure 
11,  the  most  signifieant  modifieation  is  in  regard  to  the  roughness  or  waviness  along  the 
traek’s  15.5  km  (50,800  feet).  In  reality,  what  eauses  the  sled  to  “bounee”  down  the  traek 
are  the  traek  irregularities  on  the  top  surfaee  as  well  ever  ehanging  gap  between  the  sides 
of  the  shoe  and  the  rail  in  the  width  dimension.  These  irregularities  are  mostly  a  result  of 
inherent  manufaeturing  and  rail  alignment  toleranees.  As  pointed  out  by  Cinnamon,  the 
eurrent  speeifieation  for  rail  flatness  is  no  more  than  0.0635  em  of  rail  height  ehange  in  a 
132  em  length  and  no  more  than  0.19  em  ehange  in  a  2.54  em  length  at  the  rail  seams 
[7:9-14,15].  Figure  10  illustrates  some  of  the  eonditions  that  might  exist  along  the  rail  or 
at  a  seam  yet  still  be  within  the  eurrent  HHSTT  toleranee.  Rail  eonditions  sueh  as  these, 
in  effeet,  are  part  of  what  eauses  the  sled  to  vibrate.  The  influenee  of  an  asperity  like  the 
ones  shown  below  on  wear  will  be  examined  later. 

J_ 

— ^  0.063.'  cm 

_L 

— ^  0.0635  cm 

J_ 

— ^  0.0635  cm 

_L 

0.19  cm 

Figure  10  -  Examples  of  Possible  Rail  Asperities 
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After  these  rail  irregularities  were  ineorporated,  the  model  produeed  results  that 
correlated  with  measured  dynamics  from  three  separate  sled  runs.  This  rail  geometry 
data  is  acquired  and  input  into  the  DADS  model  every  time  the  Holloman  survey  team 
physically  inspects  and  measures  the  height  and  width  of  the  rails  at  66  cm  (26  in) 
intervals  along  the  length  [10].  Asa  result  of  this  process,  the  engineers  at  Holloman 
have  great  confidence  that  DADS  accurately  predicts  the  dynamic  environment.  [12]. 


Figure  11  -  Nike  Sled  Used  in  Mission  80X-A1  [1] 


Courtesy  of  the  HHSTT,  Figure  1 1  shows  a  picture  of  Mission  80X-A1,  a  Nike 
sled  from  October  2000,  an  actual  test  that  was  part  of  the  Hypersonic  Upgrade  Program 
(HUP).  The  DADS  model  used  in  this  research  is  based  on  this  sled  configuration. 


4  3 


2 


Figure  12  -  Cross  Section  of  Rail  Showing  DADS  Contact  Locations  [1] 

DADS  reports  normal  forces  at  locations  1  through  6  in  the  rail  cross-section 
shown  in  Figure  12  for  each  shoe.  Locations  2  and  5  are  the  subject  of  future  work  for 
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the  HHSTT  to  include  the  effects  of  wind  on  yawing.  DADS  also  reports  vertical  and 
down  rail  velocities  for  the  sled  at  each  of  the  four  shoes.  This  research  will  focus  only 
on  forces  acting  on  locations  3  and  4  for  the  right  rear  shoe.  The  values  for  these  two 
locations  will  be  averaged  in  order  to  eliminate  the  rolling  motion  of  the  sled  and 
estimate  the  average  force  exerted  over  the  entire  area  of  the  shoe.  Finally,  the  particular 
DADS  file  used  in  this  research  (partially  listed  in  Appendix  D)  reports  all  of  this  data  as 
a  function  of  time  for  every  0.0001  seconds,  or  100  microseconds.  This  becomes 
important  later  in  the  subsequent  discussions  of  heat  transfer  and  wear. 


Figure  13  -  DADS  Model  for  Mission  80X-A1  Hybrid 


Mission  80X-A1*  Data  Characteristics 

Before  proceeding  further,  it  is  worth  looking  at  the  DADS  data  for  the  hybrid 
Mission  80X-A1  that  was  discussed  in  chapter  I.  This  actual  sled  only  achieved  a 
velocity  of  762  m/s  (2500  ft/s).  The  model  was  increased  for  the  purposes  of  studying 
high  velocities.  Figure  14  depicts  distance,  velocity  and  acceleration  down  the  track  as 
well  as  vertical  force  and  velocity  all  as  a  function  of  time. 
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Figure  14  -  Characteristics  of  DADS  Data 

The  distance  traveled  by  the  sled  model  increases  parabolically  with  time,  with 
the  velocity  verifying  this  by  its  linear  increase  with  time,  and  the  acceleration  is 
approximately  constant  at  1219  m/s  (4000  ft/s  ).  With  the  velocity  being  supplied  by 
DADS,  the  distance  was  found  by  integrating  the  velocity  with  respect  to  time  using  the 
trapezoid  rule.  Recall  that  the  wear  rate  spoken  about  in  equation  (2)  was  in  terms  of 
volume  of  wear  per  distance  slid;  the  distance  therefore  is  used  in  calculating  the  wear 
volume.  The  acceleration  was  found  by  (vrV/_;)/zlt.  Taking  a  closer  look  at  the  vertical 
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velocity  and  vertical  force  in  Figure  15  shows  how  the  sled  actually  comes  off  of  the  rail 
(where  the  force  goes  to  zero). 


Time  (s) 


X  10® 


Time  (s) 


Figure  15  -  DADS  Vertical  Velocity  and  Vertical  Force  for  1,515  -  1,525  s 


CTH 

Developed  by  Sandia  National  Laboratories,  CTH  is  a  hydrocode  that  was 
designed  to  simulate  large  deformation  and  shock  physics  of  multi-material,  multi¬ 
dimensional  models.  CTH  focuses  on  the  conservation  of  energy  and  momentum 
equations  along  with  equations  of  state  and  a  constitutive  model.  CTH  uses  finite 
differencing  and  finite  volume  methods  to  solve  these  equations  by  expressing  them  in 
Lagrangian  coordinates  for  the  material  reference  frames  and  Eulerian  coordinates  for  the 
spatial  frames.  For  a  more  detailed  discussion  on  CTH  solution  methods  see  Cinnamon’s 
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work  in  the  section  on  Hypervelocity  Impact  Solution  Procedure  [7:1-3].  Both  Laird  [17] 
and  Szmerekovsky  [24]  also  give  details  on  the  CTH  algorithm. 


Montgomery  and  Lim  and  Ashby 

As  previously  discussed,  the  Lim  and  Ashby  models  were  selected  because  of 
their  development  and  use  of  experimental  data.  The  Montgomery  pin-on-disk 
experiments  provide  the  sufficient  data  to  evaluate  the  Lim  and  Ashby  wear  equations.  A 
portion  of  the  experimental  pin-on-disk  data  from  Montgomery  is  in  Table  1,  while  the 
complete  data  can  be  seen  in  Appendix  E. 


Table  1  -  Excerpt  from  Montgomery  Pin-on-Disk  Experimentation  [19] 


velocity 

Bearing 

PV 

f 

fPV 

Sliding  Wi.  loss  Length 

Wear  rate 

(rts~’) 

pressure 
(lb  in"  '**) 

(X  10“®) 

(X  10' 

’*)  distance  (mg)  loss 

(ft)  (mil) 

(X  to’ 

.  3 

in  It  ) 

150 

11  300 

1.7 

0.26 

0.44 

5.7  15.0  25 

200 

600 

11  900 

7.1 

0.27 

1.92 

7.2  0.6  3 

6 

The  third  column  of  Table  1  is  the  Pv  value  discussed  earlier.  In  addition  to 
measuring  the  normal  force,  Montgomery  also  measured  the  friction  force  on  the  pins  by 
using  strain  gages  on  the  pin  holder,  which  was  also  discussed  earlier.  This  allowed  for 
the  calculation  of  the  coefficient  of  friction,  //,  from  the  equation  Ff  =  juF^.  This  value  is 
located  in  column  four  labeled  ‘/’.  The  far  right  column  is  the  experimental  wear  in 
(in  /ft)  presumably  obtained  by  measuring  the  weight  loss  and  dividing  by  the  density 
and  sliding  distance.  The  diameter  of  the  pin  used  in  these  experiments  was  0.203  cm 
(0.080  inches).  This  data  was  subjected  to  the  Lim  and  Ashby  equations  for  melt,  plastic 
(both  mild  and  severe),  and  oxidation  dominated  wear  and  plotted  along  with  the 
experimental  value  for  wear.  These  results  can  be  seen  in  Figure  16.  Note  that  oxidation 
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wear  does  not  appear  on  the  ehart  because  the  results  from  the  equation  were  not  realistic; 
they  were  extremely  large.  Additionally,  Lim  and  Ashby  did  not  use  the  Montgomery 
data  to  validate  their  oxidation-dominated  model  like  they  did  for  the  melt  wear  model. 
This  may  suggest  that  oxidation  wear  does  not  occur  at  the  speeds  considered  in 
Montgomery’s  data. 

Before  proceeding  with  discussing  the  results  of  applying  the  Lim  and  Ashby 
equations  to  the  Montgomery  data,  it  is  necessary  to  expand  on  the  a  term  used  by  Lim 
and  Ashby  in  equation  (1 1)  to  characterize  the  heat  distribution  split  between  the  pin  and 
the  disk.  This  is  important,  because  a  slightly  different  assumption  will  be  made  later 
when  looking  at  the  shoe  and  rail  for  determining  the  temperature  profile.  Lim  and 
Ashby  declared  a  to  be 

l  +  1  +  j5{m 

This  is  a  valid  expression  if 

P  =  -  (16) 

They  suggest  that  is  a  dimensionless  number  of  general  order  1”,  but  use  =  6  [19:7] 
(recalling  that  ro  is  the  pin  radius,  1.016  mm  (0.040  inches)  in  Montgomery’s  work). 

Other  important  values  are  listed  in  Table  2.  However  4,  the  “equivalent  linear  diffusion 
distance”,  is  not  so  easily  known.  If  one  should  solve  for  it  in  equation  (16),  using  /3=  6, 
lb  would  equal  6.096  mm  (0.24  inches),  and  applying  either  expression  from  equation 
(15)  into  the  melt  wear  model,  equation  (11),  yields  results  with  negative  values 
representing  growth,  not  wear.  On  the  other  hand,  if  4  =  1.016  mm  (0.04  in),  (solved  by 
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using  j5=\,  but  still  using  6  in  equation  (11))  and  the  left  expression  of  equation  (15) 
is  used  in  equation  (1 1)  to  obtain  the  melt  wear  values,  the  melt  wear  results  in  what  is 
shown  in  Figure  16. 

Table  2  -  Values  used  for  Lim  and  Ashby  Equations  for  the  Pin-on-Disk  (Steel) 


Value 

Description 

SI 

English 

ro 

Pin  Radius 

1.016  mm 

0.040  in 

H 

Pin  Flardness 

1.00  GPa 

1.45  X  10^  psi 

a 

Thermal  Diffusivity 

9.1  X  lO'Ws 

0.014  in^/s 

K 

Thermal  Conduetivity 

0.41  W/emK 

5.12  Ibf/sR 

To 

Pin  Temperature 

297  K 

535  R 

Tm 

Melt  Temperature 

1800  K 

3240  R 

T 

Thermal  Quantity 

222  K 

400  R 

L 

Latent  Pleat  of  Melting 

2.1  X  10“^  jW 

3.05  X  10^  psi 

r*  =  aH/K 
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Figure  16  -  Montgomery  Experimental  Data  and  Montgomery  Data  used  in  Lim  and 

Ashby  Equations 

•2 

The  wear  rate  is  plotted  as  volume  (mm  )  per  sliding  distance  (m)  versus  the 
quantity  of  pressure  multiplied  by  the  velocity  (GPa  m/s).  Two  conclusions  can  be  drawn 
from  Figure  16;  first,  melt  wear  is  present  if  the  correct  expression  for  a  is  used  in 


equation  (11),  and  the  correct  value  of  P  and  4  are  chosen.  And  second,  the  plastic 
(mechanical)  wear  model,  as  defined  earlier,  can  be  a  viable  solution  in  defining  wear 
which  is  simply  a  function  of  the  non-dimensional  normal  force  multiplied  by  the 
constant  kA,  as  shown  in  equation  (12). 
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Further  developing  the  second  conclusion,  kA  being  somewhat  arbitrary,  can  be 
solved  for  using  the  empirical  data  given  by  Montgomery  in  the  following  manner  to 
provide  a  mechanical  wear  function  for  all  values  of  Pv  within  the  experimental  range. 


= 


w 


Montgomery 


H 


A  P 

Montgomery 


(17) 


Here,  WMontgomery  and  P Montgomery  arc  thc  experimental  wear  rate  and  pressure  respectively 
from  Montgomery.  These  values  of  kA  are  shown  in  Figure  17  as  ''kA  \  Notice  that  some 
values  of  kA  at  low  Pv  values  are  out  of  the  suggested  range  of  5  x  10'  to  5  x  10'  given 
by  Lim  and  Ashby  in  the  section  on  plastic  (mechanical)  wear.  These  values  were 
manually  adjusted  to  the  range  limit  of  5  x  10'  suggested  by  Lim  and  Ashby.  These 
numbers  are  represented  in  Figure  17  as  adjusted”.  Finally,  these  adjusted  values 
were  averaged  over  three  arbitrary  ranges  of  Pv,  which  can  be  seen  in  Figure  17  as  the 
series  averaged”. 
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Figure  17  -  Adjusted  Ra  Values  for  Montgomery  Data 

This  final  set  of  kA  values  for  the  three  different  Pv  ranges,  shown  in  Table  3, 
were  then  used  to  calculate  the  theoretical  wear  as  defined  by  equation  (4).  Figure  18 
shows  how  this  method  correlates  with  the  experimental  data. 

Table  3  -  Values  for  Ra  Across  Pv  Range 


Pv 

Ua 

(GPa  m/s) 

(-) 

<4.2 

3.53  X  10'^ 

4.2-31.5 

8.16  X  10'^ 

>31.5 

1.30  X  10'^ 

35 


5 


4.5 
4 

3.5 

"e 

ct'  3 

E 

<u  2.5 

03 

a: 

S-  O 
03  ^ 

0 

§ 

1.5 
1 

0.5 
0 

0  10  20  30  40  50  60  70 

Pv(GPa  m/s) 

Figure  18  -  Montgomery  Experimental  Data  with  Calculated  Wear 
Using  Equation  (17)  and  Adjusted  kA  Values  from  Table  3 

The  method  for  finding  the  appropriate  kA  eould  be  modified  to  use  a  least  squares 
fit  or  some  eurve  fit  function  instead  of  using  an  average  in  three  different  arbitrary 
regions.  For  the  current  effort,  the  values  in  Table  3  will  be  used  to  calculate  wear  for  the 
Holloman  sled.  This  leads  to  the  next  step,  the  application  of  these  wear  models, 
equation  (11)  and  equation  (12)  to  the  DADS  data  to  calculate  wear. 

DADS  and  Lim  and  Ashby 

Just  as  Table  2  defined  some  necessary  values  for  steel  in  the  pin-on-disk 
experimentation.  Table  4  provides  the  same  data  for  VascoMax  as  related  to  the  HHSTT. 
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Table  4  -  Values  used  for  Lim  and  Ashby  Equations  for  the  HHSTT  (VascoMax) 


Value 

Description 

SI 

English 

ro 

Equivalent  Radius 

50.80  mm 

2.0  in 

Aft 

Shoe  Area 

206.45  cm^ 

32.0  in^ 

Hv 

Shoe  Hardness 

1.00  GPa 

1.45  X  10^  psi 

Qy 

Thermal  Diffusivity 

4.5  X  10'^  m^/s 

0.007  in^/s 

Kv 

Thermal  Conductivity 

0.31  W/cmK 

3.85  Ibf/sR 

CpV 

Specific  Heat 

0.858  J/g  K 

7.39  X  10^  inWR 

To 

Shoe  Temperature 

297  K 

535  R 

TmV 

Melt  Temperature 

1685  K 

3033  R 

Tv 

Thermal  Quantity 

146  K 

264  R 

Lv 

Latent  Heat  for  Melting 

2.1  X  10‘^  jW 

3.05  X  10^  psi 

Subscript  V  denotes  VascoMax 

The  thermal  eonduetivity,  thermal  diffusivity  and  specific  heat  for  Vascomax  were  taken 
from  Laird’s  work  where  he  assumed  constant  properties  at  a  temperature  of  700  K 
[17:84],  The  material  hardness  and  latent  heat  of  melting  are  the  same  as  that  of  steel. 

The  previous  section  concluded  that  equation  (11),  with  a  properly  chosen  value 
for  lb  and  equation  (12),  with  the  appropriate  value,  would  provide  analytical  results 
for  predicting  melt  and  mechanical  wear  respectively  for  the  shoe.  In  order  to  use 
equation  (11),  however,  the  coefficient  of  friction  is  needed. 
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Coefficient  of  Friction 


Montgomery  measured  the  eoeffieient  of  friction,  ju.  DADS,  unfortunately,  does 
not  provide  this  information.  Lim  and  Ashby  used  the  following  equation  for  their 
models  which  is  based  on  a  variety  of  experimental  data: 

//  =  0.78-0. 13  logio(v)  (18) 

/u  is  expressed  solely  as  a  function  of  normalized  velocity  [18:7].  However, 
when  V  reaches  10^,  //  equals  zero,  and  any  further  increase  in  v  causes  //  to  become 
negative.  For  the  Holloman  sled,  v  reaches  10^  or  88.9  m/s  (291.6  ft/s)  in  the  first  0.07 
seconds  of  the  run  according  to  equation  (6)  and  Figure  14.  Therefore,  the  coefficient  of 
friction,  //,  was  determined  by  using  a  curve  fit  from  the  Montgomery  pin-on-disk  data, 

not  only  as  a  function  of  v  but  P  {ox  F  )  as  well.  Pv  is  used  instead  of  Pv  so  that  it  can 
be  used  in  a  power  function  without  having  to  deal  with  the  units.  Generating  a  curve  fit 
of  the  plotted  coefficient  of  friction  vs.  normalized  Pv  shows  a  power  function  that  can 
also  be  applied  to  the  DADS  data  for  the  range  of  normalized  Pv  under  examination. 

Figure  19  shows  the  Montgomery  friction  data  plotted  against  Pv  along  with  two  curves, 
the  higher  curve  is  representative  of  the  function  shown,  while  the  lower  curve  has  been 
modified  so  that  the  coefficient  of  friction  will  asymptotically  approach  0.02  for  the  range 
of  normalized  Pv  shown  (0  -  8000).  In  a  separate  paper,  Montgomery  recognized  that 
experimental  pin-on-disk  data  for  the  coefficient  of  friction  was  an  order  of  magnitude 
higher  than  the  actual  value.  His  work  suggested  that  the  true  value  of  the  coefficient  of 
friction  approaches  0.02  as  Pv  increases  [20:242]. 
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Figure  19  -  Function  for  Coefficient  of  Friction  vs  Pv  Based  on  Montgomery  Data 

The  modified  power  function  was  simply  determined  by  subtracting  a  constant 
from  the  actual  curve.  The  constant  was  calculated  by  taking  the  minimum  value  (for  the 

Pv  range  of  0  -  8000)  on  the  actual  curve  and  subtracting  0.02. 

//  =  1.5704(Pv)-'’2299 

where  (19) 

5  =  min{l.5704(Pv)  }- 0.02 

Equation  (19)  was  then  used  to  calculate  the  coefficient  of  friction  for  the  DADS  data. 
The  results  are  displayed  in  Figure  20.  These  values  will  be  used  in  subsequent 
calculations  for  melt  wear  as  well  as  frictional  heat  generation. 
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Figure  20  -  Coefficient  of  Friction  vs.  Normalized  Pv  for  DADS  Data 

Calculation  of  Wear  for  DADS 

In  order  to  calculate  wear  from  equation  (12)  for  the  DADS  data,F  was  calculated 
from  equation  (5),  and  Pv  was  calculated  by  dividing  the  normal  force  by  the  nominal 
area  of  the  shoe,  An,  assuming  uniform  contact,  and  multiplying  by  the  velocity. 

Knowing  the  Pv  term,  Ica  for  each  time  step  can  be  taken  from  Table  3  and  be  used  to 
calculate  the  wear  for  that  time  step.  Similarly,  equation  (11)  will  provide  the  melt  wear 
quantities,  but  first  a  needs  to  be  calculated  from  the  left-hand  expression  in  equation 
(15).  The  value  of  4  was  chosen  as  1.016  cm  (0.04  in),  and  f=  6,  similar  to  the  Lim  and 
Ashby  model  that  produced  the  melt  wear  results  shown  in  Figure  16. 
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This  leads  to  two  quantities,  wear  due  to  melting  from  equation  (11)  and  wear  as  a 
function  of  kA  from  equation  (12),  both  as  a  function  of  time,  or  for  that  matter  a  function 
of  Pv  or  distance.  For  the  remainder  of  this  research,  these  quantities  will  be  referred  to 
as  (mechanical)  and  Wmeit-  Figure  21  shows  and  Wmeit  as  a  function  of  distance. 

The  second  plot  of  Figure  21  is  a  close-up  look  over  a  distance  of  50  m  centered  at  the 
600  m  mark.  This  illustrates  the  “bouncing”  effect  of  the  sled.  When  the  wear  rates  goes 
to  zero,  the  sled  is  not  in  contact  with  the  rail.  By  taking  the  area  under  these  curves 
(using  the  trapezoid  rule  for  integration),  it  is  possible  to  determine  the  total  wear  volume 
and  a  wear  thickness  in  the  shoe  due  to  each  mechanism.  These  results  will  be  discussed 
in  Chapter  III. 
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Figure  21  -  Wear  Rates  for  DADS  Data  for  Entire  Sled  Run  and  for  575  -  625  m 


Heat  Transfer 

The  DADS  supplied  normal  force  is  used  to  calculate  the  amount  of  heat  transfer 
that  takes  place  in  the  shoe  due  to  the  frictional  heating  at  the  boundary  of  the  shoe  and 
rail  interface.  This  method  was  first  performed  by  Laird  in  his  work  on  gouging  [17:83- 
89].  For  the  purposes  of  this  research,  a  one-dimensional  heat  transfer  analysis  was 
considered  sufficient.  Later  research  should  consider  the  two-dimensional  effects  as  well 
as  other  heat  sources  such  as  the  stagnation  temperatures  at  the  front  of  the  shoe  resulting 
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from  the  hypersonic  flow.  Frictional  heating  can  be  expressed  as  a  function  of  time,  in 
terms  of  power  by  the  following  equation; 

q{t)  =  ^{t)F ^{t)v{t)  (20) 

The  assumption  is  made  that  because  the  thermal  conductivities  of  the  1080  steel 
rail  and  the  VacoMax  300  shoe,  an  alloy  of  steel,  are  very  similar,  50%  of  the  heat 
generated  by  friction  is  transferred  to  the  rail  while  the  remaining  50%  is  transferred  to 
the  shoe.  Finally  the  heat  transfer  into  the  shoe  in  terms  of  flux  (W/m  )  can  be 
determined  for  each  time  step  of  the  DADS  data  by  using  the  following  expression: 

(21) 

n 

The  %  represents  the  heat  transfer  split  between  the  rail  and  shoe.  Laird,  in  his 
work,  solved  a  heat  equation  with  a  constant  temperature  boundary  condition  at  the 
rail/shoe  interface  determined  from  the  total  energy  generated  by  frictional  heating  over 
the  entire  run  using  the  average  normal  force  value. 
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Figure  22  -  Flux  and  Temperature  vs  Time  for  the  First  0,3  Seconds 

In  this  case,  the  flux,  q  ”(t),  will  be  used  as  the  boundary  condition  for  each  time  step  in 
the  solution.  The  flux  as  a  function  of  time,  for  the  first  0.3  seconds  of  the  run  is 
displayed  in  Figure  22  along  with  the  corresponding  temperature. 

The  temperature  through  the  thickness  of  the  shoe  as  a  function  of  time  is 
determined  by  solving  the  heat  equation  in  one  dimension; 

dT 


dt 


=  a. 


( 


t>0,  y>  0 


(22) 


where  y  is  the  positive  direction  into  the  shoe.  The  one-dimensional  solution  was  thought 
to  be  sufficient  for  determining  the  temperature  profile  in  the  shoe  especially  for  the  sled 
traveling  at  slower  speeds.  A  more  accurate,  two-dimensional  approach  should  be 
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investigated  at  higher  sled  speeds.  The  a ^  term  in  equation  (22)  is  the  thermal  diffusivity 
defined  by 


Uy 


K 

pc 


V 


Pv 


(23) 


Ky  is  the  thermal  eonductivity,  p  is  the  density  and  Cpv  is  the  speeific  heat  of  VaseoMax. 

A  finite  differeneing  method  was  used  to  solve  this  one-dimensional  heat  equation 
explieitly  using  a  forward  differenee  on  time  and  eentral  difference  on  space  by 


T"'  =  T-  +  a,  (t::,  --it- +  T-_, )  (24) 

At 

i  represents  the  spatial  coordinate,  and  n  represents  time.  If  a  second  order  approximation 
is  used  for  the  boundary,  the  flux  condition 


dT  _  q'\t) 

dy  Ky 

can  be  represented  for  each  time  step  as 


(25) 


7;  =  7;  +  2Ay 


q\t) 

K,, 


(26) 


Where  To  is  at  an  imaginary  point  outside  the  boundary  and  is  used  in  equation  (24)  in  the 
Tyi  position  to  calculate  the  temperature  at  the  boundary.  Figure  23  illustrates  the  flux 
boundary  condition.  This  indicates  that  the  slope,  dT/dy,  at  the  boundary  is  negative 
according  to  the  right  side  of  equation  (25),  and  therefore  the  temperature,  moving  away 
from  the  boundary,  is  decreasing. 
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Figure  23  -  Heat  Diffusion  with  a  Flux  Boundary  Condition  [15:63] 

For  this  explicit  method,  it  is  important  to  check  the  stability  of  the  solution  by 
maintaining  that  the  quantity  ,  is  less  than  0.5  [6:702],  Because  the  heat  equation 

is  second  order  in  the  spatial  dimension  and  first  order  in  time,  the  solution  requires  a 
second  boundary  condition,  as  well  as  an  initial  condition.  Therefore  the  boundary 
condition  in  the  far  field  was  a  shoe  temperature  equal  to  ambient,  or  297  K.  The  initial 
condition  was  a  uniform  ambient  temperature  throughout  the  shoe. 

The  CTH  Plane  Strain  Model 

The  numerical  model  being  used  for  high  velocity  wear  was  developed  by  Laird 
[17]  and  improved  upon  by  both  Szmerekovsky  [24]  and  later  Cinnamon  [7]. 
Szmerekovsky  laid  the  groundwork  for  the  two  dimensional  plane-strain  model  by 
conducting  a  mesh  convergence  study  and  computing  the  simulation  run-time  limit.  The 
mesh  converged  at  a  square  cell  dimension  of  0.002  cm.  Using  a  smaller  mesh  would 
increase  total  simulation  time  without  significantly  improving  the  results.  The  reflected 
stress  wave  at  the  material  interface  boundary  limits  the  simulation  to  less  than  20 
microseconds.  A  10  microsecond  simulation  time  is  used  to  ensure  there  would  be  no 


46 


interference  from  the  wave  reflection.  Cinnamon  added  the  strain  rate  dependent  Zerilli- 
Armstrong  constitutive  model  for  VascoMax  300  and  1080  steel  [7:3-1 1  to  3-12]. 

a  =  A  +  (c^  +  Cj  +  c^s:"  (27) 


The  plane  strain  model  developed  by  Szmerekovsky  is  described  in  Figure  24. 
The  intention  is  to  model  a  real  world,  three-dimensional  scenario  in  a  two  dimensional 
plane  strain  simulation,  taking  care  to  maintain  the  kinetic  energy  and  momentum.  The 
model  represents  a  sled  of  841  kg  (1854  Ibm).  The  mass  is  evenly  distributed  over  the 
four  shoes  and  sliced  from  10.16  cm  (4  in)  to  1  cm  (0.394  in)  to  make  it  a  plane  strain 
model  (CTH  uses  length  units  of  centimeters).  Platinum,  being  the  most  dense  material 
defined  in  CTH,  /?  =  21 .45  gm/cm  ,  is  used  for  the  sled  mass  having  a  total  height  of 


46.53  cm  (18.32  inches).  This  calculation  is  detailed  in  Appendix  C. 
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Figure  24  -  Plane  Strain  CTH  Model  Development  [24:5-2] 


47 


It  is  important  to  note  that  this  method  neglects  the  potential  effects  of  sled  roll, 
pitch  and  yaw,  and  the  possibility  of  the  weight  being  unevenly  distributed  on  the  sled. 
These  considerations  should  be  researched  in  future  efforts. 

Correlating  CTH  with  DADS 

In  order  to  use  CTH  to  model  what  DADS  predicts,  the  attempt  was  made  to 
associate  the  DADS  normal  force  with  the  forces  experienced  in  the  CTH  model  at  a 
particular  horizontal  velocity.  Holloman  historically  suggested  using  a  vertical  velocity 
of-1  m/s  for  all  cases,  because  this  is  roughly  the  average  vertical  velocity  produced  by 
DADS  over  an  entire  sled  run.  Previous  research  involving  gouging  has  been  run  with 
this  -1  m/s  vertical  velocity.  Because  CTH  does  not  output  force,  an  attempt  was  made 
to  determine  the  downward  force  in  CTH  by  looking  at  the  rigid  body  motion  over  a 
chosen  time  interval.  Initially  this  was  performed  at  100  m/s  sled  velocity.  The  rigid 
body  was  identified  as  the  concentration  of  the  sled  mass  that  had  the  same  vertical 
velocity  from  one  time  period  to  the  next,  suggesting  no  deformation.  The  force  was 
determined  with  Newton’s  second  law  by  estimating  the  mass  of  the  rigid  body  from  the 
dimensions  in  Figure  25  to  obtain  the  volume,  multiplying  by  density  and  then 
multiplying  by  the  change  in  the  vertical  velocity  over  time. 
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Figure  25  -  CTH  Vertical  Force  Prediction  Using  Av/At 

Tracer  point  39,  which  lies  within  the  rigid  body,  hinted  toward  the  fact  that  the  rigid 
body  was  accelerating  downward  as  indicated  by  the  chart  on  the  right.  The  scale  of  the 
chart  on  the  left  provided  a  coarse  estimation  of  the  velocity.  This  calculation  resulted  in 
an  acceleration  of  -21.33  m/s  and  a  force  of  3,123  N,  remarkably  close  to  the  3,363  N 
that  DADS  supplied  at  this  speed  for  this  particular  sled  file.  This  technique  was  also 
applied  to  the  sled  traveling  at  1000  m/s.  The  results  were  nearly  identical  to  the  100  m/s 
case,  suggesting  that  this  technique  was  inadequate,  because  DADS  produced  a  much 
larger  vertical  force  for  the  higher  speed.  Consequently,  the  acceleration  in  CTH  that  the 
sled  experienced  was  mostly  a  product  of  the  acceleration  due  to  gravity.  What  resulted 
from  this  effort  is  that  the  gravity  term  in  CTH  was  discovered.  Previous  models  did  not 
take  gravity  into  consideration,  however  previous  models  did  not  simulate  the  sled  at  the 
relatively  low  speed  of  100  m/s  where  gravity  is  more  of  a  factor  because  the  system’s 
momentum  is  lower.  The  21 .33  m/s  acceleration  in  retrospect  was  well  within  the  error 
tolerance  of  this  calculation  to  sufficiently  say  it  was  in  fact  gravity.  The  conclusion  was 
that  CTH  focuses  more  on  momentum  conservation  than  forces  in  the  sense  that  forces 
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are  neither  input  nor  output  parameters  whereas  veloeities  are  both  input  and  output 
parameters.  The  author  reeognizes  that  foree  and  veloeity  are  related  through  the  ehange 
in  momentum  over  time.  Therefore  the  eorrelation  between  CTH  and  DADS  would  be 
suffieiently  satisfied  if,  instead  of  using  -1  m/s  vertieal  veloeity  the  aetual  DADS 
produeed  vertieal  veloeity  was  used.  Reeall  from  Figure  14,  that  DADS  ealeulates  a 
vertieal  veloeity  as  well  as  a  horizontal  sled  veloeity  as  a  funetion  of  time.  The  teehnique 
would  involve  ehoosing  the  simulation  veloeity  then  determining  the  eorresponding 
vertical  velocity. 

Defining  the  Mesh 

As  mentioned  in  the  section  on  the  plane  strain  model,  the  smallest  mesh  size  is 
0.002  cm.  Considering  that  computational  time  is  a  function  of  the  number  of  cells  and 
therefore  the  cell  size,  it  is  more  efficient  to  only  use  the  0.002  cm  cell  size  in  areas  of 
interest,  namely  at  the  material  interfaces.  It  is  also  important  to  maintain  an  aspect  ratio 
close  to  unity  for  all  cells.  The  mesh  definition  that  was  used  for  the  final  analysis  can  be 
seen  in  Appendix  B,  CTH  Input  File,  while  a  typical  model  can  be  seen  in  Figure  26. 

The  figure  on  the  left  is  depicting  the  mesh.  The  figure  on  the  right  is  material  definition 
at  time  equals  zero.  Overlaying  the  two  figures  gives  an  indication  of  the  where  the 
model  fits  in  the  mesh.  The  idea  was  to  use  a  fine  mesh  in  the  places  of  interest,  namely 
the  shoe/rail  interface  and  the  leading  edge  of  the  shoe.  The  typical  runtime  for  a 
simulation  using  this  mesh  was  roughly  eight  hours. 
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Figure  26  -  Representative  Mesh 


CTH  Material  Interface 

Cinnamon  conducted  a  study  to  determine  which  of  the  three  material  interface 
scenarios  within  CTH  would  be  best  for  the  gouging  problem.  The  three  techniques  used 
are  no-slide,  slide-line  and  boundary  layer  approaches.  The  slide-line  technique,  which 
sets  the  shear  strengths  to  zero  to  allow  sliding  to  occur,  generated  instabilities  within  the 
solution  in  that  shear  stresses  developed  far  away  from  the  interface.  The  boundary  layer 
method  was  found  to  be  valid  only  for  two-dimensional  axi-symmetric  problems  using 
one  computational  processor  to  solve  the  problem  compiled  by  CTH.  Therefore  the  no¬ 
slide  method  was  chosen  as  the  best  method  to  model  the  material  interface  for  the 
gouging  problem  [7:6-30].  The  same  method  was  used  for  the  wear  problem. 

The  concept  of  no-slip  suggests  that  no  sliding  will  occur  during  the  simulation. 
However,  CTH  handles  the  problem  of  sliding  under  no-slip  by  requiring  the  materials  to 
fail  under  shear.  Shearing  is  brought  about  by  giving  the  sled  and  rail  different  velocities. 
Initially  all  materials  are  considered  joined  at  the  beginning  of  each  time  step  where 
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contact  exists.  If  the  ease  oeeurs  where  a  cell  would  have  a  mixture  of  one  or  more 
materials,  the  cell  strength  is  governed  by  the  material  volume  fraetions.  As  it  turns  out, 
this  is  a  good  representation  of  frietion  on  the  miero-struetural  scale  discussed  in  Chapter 
I.  This  motion  can  generate  plasticity,  heat,  material  fracture  and  debris,  all  of  whieh 
eontribute  to  wear  direetly  or  indirectly.  In  addition,  the  use  of  a  lower  strength,  thin 
layered  coating  on  the  rail,  sueh  as  epoxy,  gives  the  effective  sliding  boundary  condition 
as  will  be  discussed  along  with  the  result  in  the  next  seetion. 


52 


Ill,  Results  and  Analysis 


Melt  Wear 

Lim  and  Ashby  state  that  “Agreement  is  good”  between  the  experimental  data 
(they  were  referring  to  the  Montgomery  data)  and  their  wear  model.  [18:13].  Figure  27 
is  the  graphieal  representation  of  their  wear  model  along  with  the  experimental  data 
plotted  in  3D  spaee  {v,  F  ,w).  The  contours  are  the  lines  established  by  the  model 
(equation  (11))  and  represent  lines  of  constant  wear  rate.  The  numbers  next  to  the 
experimental  data  points  are  logio(w). 

SLIDING  VELOCITY  v  (m/s) 


NORMALISED  VELOCITY  v 

Figure  27  -  Lim  and  Ashby  Model  for  Melt  Wear  [18:13] 


53 


However,  a  eloser  look  in  three  dimensions  shows  that  the  values  for  iv  are  elose  but  not 
within  the  contours  (Figure  28).  This  would  indicate  that  the  model  is  not  entirely 
accurate. 


Figure  28  -  Three  Dimensional  View  of  Lim  and  Ashby’s  Melt  Wear  Model 

Accordingly,  when  the  melt  wear  model  was  applied  to  the  DADS  data  as  discussed 
previously,  acceptable  values  for  melt  wear  were  achieved  using  the  same  values  for  (3 
and  h  that  were  used  for  the  Montgomery  data,  that  is  6  and  4  =  1.016  mm  (0.04  in). 
The  results  shown  in  and  Figure  30  are  similar,  however  the  DADS  data  points  are  not  as 
close  to  the  contour  lines  as  the  Montgomery  data  was.  This  is  because  the  function  used 
in  the  model  for  the  coefficient  of  friction,  //,  is  different  than  the  function  used  for  the 
calculation  of  melt  wear  using  the  DADS  data  (equation  (18)  versus  equation  (19)).  Also 
recall  that  the  geometry  used  to  create  the  model  was  based  on  the  pin-on-disk  scenario 
with  a  very  small  pin  diameter  whereas  the  DADS  data  uses  a  much  larger  radius  to 
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represent  the  shoe.  However  this  could  simply  be  an  indication  that  either  the  model  is 
invalid  or  the  shoe  does  not  experience  melt  wear. 


Figure  29  -  Melt  Wear  Results  for  DADS  Data 
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Figure  30  -  Three  Dimensional  View  of  Melt  Wear  for  DADS  Data 


Lim  and  Ashby  Equations  Predicting  Wear 


A  MATLAB  code  was  constructed  to  model  the  shoe  wear  from  the  DADS  data 


using  the  Lim  and  Ashby  equation  for  wear  wm  and  Wme/t-  This  was  detailed  in  the  last 
chapter.  For  the  numerical  modeling  in  CTH,  two  veloeities  were  ehosen  to  be  analyzed 
in  detail  with  one  additional  velocity  chosen  to  get  a  data  point  at  an  even  higher  velocity. 
The  velocities  were  somewhat  arbitrarily  seleeted  based  on  their  relatively  low  magnitude 
and  corresponding  vertical  velocity.  All  of  the  work  that  was  eonducted  by  Cinnamon 
with  regard  to  gouging  was  performed  at  3,000  m/s  [7].  This  researeh  focuses  more  on 
how  velocity  affects  the  shoe  by  examining  the  behavior  at  394  m/s  and  798  m/s,  in 
combination  with  several  other  parameters  and  one  case  with  a  velocity  of  1,41 1  m/s. 
Table  5  and  Figure  31  details  the  relevant  information  about  the  points  chosen  from  the 
DADS  data  to  be  analyzed. 
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Figure  31  -  Plots  Indicating  Velocities  Used  in  CTH 

The  next  sections,  including  Table  6,  discuss  the  results  from  the  analytical 
models  with  the  assumption,  as  stated  earlier,  that  the  shoe  makes  uniform  contact  with 
the  rail  over  the  entire  206.5  cm  (32  in  )  area.  This  assumption  allows  for  a  simple 
calculation  of  wear  (or  material  removal)  by  dividing  the  volume  of  wear  by  the  area. 
Recall  that  the  volume  of  wear  was  determined  by  integrating  the  wear  rate  (mm  /m)  with 
respect  to  the  distance  the  sled  traveled. 


394  m/s 

The  MATLAB  code  was  designed  to  calculate  the  wear  up  to  a  specified  input 
velocity.  The  amount  of  mechanical  wear  generated  by  the  sled  according  to  equation 
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(12)  up  to  394  m/s  was  0.0035  cm,  and  the  amount  of  melt  wear  was  0.0047  em  after 
traveling  a  distanee  of  61.1  m. 

The  temperature  distribution  resulting  from  the  frietional  heating  that  was  used  to 
put  into  CTH  is  shown  in  Figure  32.  Reeall  that  the  flux,  q”(t),  is  a  funetion  of  time,  and 
the  eoeffieient  of  Motion  is  a  funetion  of  foroe  and  velooity  as  stated  by  equation  (19). 
Depending  on  the  length  and  duration  of  the  oontaot,  a  speoifio  quantity  of  heat  is 
generated  due  to  friotion  that  aooumulates  over  time.  The  figure  on  the  left  therefore,  is 
the  ourve  that  MATLAB  generated  for  the  temperature  distribution  in  the  shoe  at  the 
point  when  the  sled  reaohes  394  m/s.  The  indioated  points  were  taken  and  used  to  oreate 
a  temperature  profile  in  CTH.  The  figure  on  the  right  is  the  representative  profile  in 
CTH.  The  maximum  temperature  seen  by  the  shoe  was  oaloulated  to  be  567  K,  whioh 
may  seem  oontrary  to  Figure  32.  Remembering  that  the  flux  varies  with  time  and  at  some 
instanoes  equals  zero,  suggests  that  the  time  the  sled  reaehes  394  m/s  or  0.31  seeonds  into 
the  run  the  sled  is  not  in  oontaot  with  the  rail,  and  heat  is  not  being  transferred  into  the 
shoe. 
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Figure  32  -  Temperature  Distribution  for  Heat  Generated  from  Friction  for  the  Sled 
up  to  394  m/s  -  MATLAB  and  CTH  Curves 


798  m/s 

Identical  to  the  previous  case,  the  MATLAB  code  was  executed  for  the  sled’s 
travel  up  to  798  m/s.  The  resulting  reduction  in  thickness  of  the  shoe  due  to  mechanical 
wear  was  0.0156  cm  and  0.01 13  cm  for  melt  wear.  The  maximum  temperature  reached 
was  854  K  in  a  distance  of  252.0  m.  The  temperature  distribution  from  frictional  heating 
for  this  speed  is  in  Figure  33.  The  MATLAB  generated  profile  is  on  the  left  along  with 
the  points  used  in  CTH.  The  CTH  profile  from  those  points  is  shown  in  the  figure  on  the 
right. 
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Figure  33  -  Temperature  Distribution  for  Heat  Generated  from  Friction  for  the  Sled 
up  to  798  m/s  -  MATLAB  and  CTH  Curves 


1,411  m/s 

The  analysis  at  1,41 1  m/s  was  not  conducted  in  as  much  detail  as  the  slower 
velocities.  The  analytically  derived  melt  wear  amount  of  0.0246  cm  and  a  mechanical 
wear  amount  of  0.0576  cm  was  calculated  for  a  distance  traveled  of  796.5  m. 

One  simulation  was  conducted  in  CTH  at  this  velocity  against  an  asperity,  at 
ambient  initial  temperature  conditions  and  without  a  coating.  For  this  reason  a 
temperature  profile  was  not  established  to  be  input  into  CTH.  However,  according  the 
heat  transfer  portion  of  the  MATLAB  code,  the  maximum  temperature  was  1451  K. 

Table  6  is  a  complete  summary  of  the  results  for  both  melt  wear  and  mechanical 
wear  calculated  from  the  analytical  model  developed  in  MATLAB.  This  code  can  be 
seen  in  Appendix  A.  Notice  that  the  last  row  of  data  is  for  the  complete  run.  This  shows 
a  total  wear  thickness  of  0.273 1  cm  for  mechanical  wear  and  0.0797  cm  for  melt  wear. 
These  values,  if  superimposed,  amount  to  approximately  0.35  cm  of  total  wear.  Although 
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these  results  may  not  be  eompletely  accurate,  they  are  realistic  considering  the  shoe  is 
2.54  cm  thick. 

Table  6  -  Analytical  Model  Results  for  Melt  and  Mechanical  Wear  and  Maximum 


Temperature 


1'^^ . 

j'M  *  ^ 

‘1  j  ^ 

r  Maximum 

Temperature 

1 1 

Melt 

Wear 

(cm) 

Mechanical 

Wear 

(cm) 

394  m/s 

570 

0.0047 

0.0035 

798  m/s 

854 

0.0113 

0.0156 

1,411  m/s 

1,451 

0.0246 

0.0576 

3,030  m/s 

3,925 

0.0797 

0.2731 

In  addition,  Figure  34  shows  how  the  DADS  data  looks  in  terms  of  wear  rate  vs. 
Pv.  Comparing  this  figure  with  Figure  16  from  the  Montgomery  data  shows  that  the  wear 
rates  are  significantly  higher  for  the  DADS  data,  however  all  of  the  DADS  data  occurs  at 
lower  Pv  values  compared  to  the  Montgomery  data. 
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CTH 

In  the  effort  to  numerieally  generate  wear  results  using  CTH,  several  scenarios  of 
the  sled/rail  system  were  simulated.  Four  different  parameters  with  two  possible 
variations  each,  were  simulated  in  all  possible  combinations  resulting  in  16  cases.  The 
parameters  are  velocity,  394  m/s  and  798  m/s  as  discussed  above;  a  flat  rail  or  a  rail  with 
an  asperity;  a  rail  with  or  without  a  coating;  and  the  rail  at  ambient  temperature  and  the 
rail  with  the  profile  generated  from  the  MATLAB  code.  A  0.01  cm  coating  was 
customarily  used  in  CTH  gouging  simulations  to  more  closely  approximate  the  sliding 
phenomenon  that  CTH  does  not  adequately  model  as  previously  discussed.  The  use  of  a 
coating  provided  excellent  results  correlation  for  past  researchers,  therefore  it  has  been 
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carried  forward  for  the  modeling  of  wear.  In  reality,  a  thin  epoxy  layer  is  frequently  used 
at  Holloman,  therefore  using  a  coating  in  the  model  is  not  impracticable. 

It  was  antieipated  that  eharacterizing  wear  in  CTH  would  be  possible  in  the 
meehanisms  of  melt  and  loealized  fracture  or  failure.  With  this  in  mind,  the  CTH  results 
were  analyzed  by  examining  the  shoe  temperature  and  the  maximum  prineipal  stress  in 
the  shoe.  The  temperature  field  is  generated  during  the  CTH  model  exeeution  from 
plastie  deformation  in  the  materials;  this  is  also  referred  to  as  plastieity.  Melt  was 
determined  by  comparing  the  temperature  field  produced  by  CTH  to  the  melt  temperature 
of  VascoMax  300  which  is  1685  K.  Localized  failure  occurs  when  the  maximum 
principal  stress  exeeeds  the  fraeture  strength  defined  for  VascoMax  300.  The  maximum 
principle  stress  is  the  largest  normal  stress  in  a  prineipal  direction  for  any  given  state  of 
stress.  The  principal  directions,  or  principal  planes,  are  those  for  which  a  given  state  of 
stress,  including  normal  and  shearing  stresses,  are  transformed  to  only  normal  stresses. 
Cinnamon  experimentally  determined  the  fraeture  strength  for  VascoMax  300  to  be  2.5 
GPa  (2.5  X  10  dyne/cm  )  [7:5-13].  CTH  eonveniently  supplies  the  maximum  principal 
stress  as  an  output  parameter. 

The  following  sections  will  examine  the  output  from  CTH  for  the  16  cases 
simulated.  The  proeedure  was  to  first  visually  find  the  highest  temperature  or  maximum 
principal  stress  from  the  respective  contours  and  then  produce  a  one-dimensional  plot  of 
that  parameter  versus  the  y  dimension,  or  through  the  thickness  of  the  shoe,  as  it  passes 
through  that  maximum  area  identified  on  the  eontour. 


63 


Flat  Rail,  394  m/s,  Coating 


For  this  case,  neither  the  maximum  principal  stress  nor  the  temperature  exceeded 
the  fracture  strength  or  the  melt  temperature  respectively.  The  maximum  temperature 
reached  in  the  shoe  material  was  approximately  329  K  at  1  microsecond,  near  the  back  of 
the  shoe,  well  short  of  melt.  Figure  35  shows  a  2D  temperature  profile  for  the  rail  shoe 
interface  as  well  as  a  ID  look  at  the  temperature  into  the  shoe  at  an  x  position  of  2.38  cm, 
both  of  which  are  for  simulation  time  equal  to  1  microsecond. 
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Figure  35  -  Temperature  for  394  m/s,  Coating,  Initial  Ambient  Temperature 


Similarly,  Figure  36  shows  areas  of  high  stress  near  the  front  of  the  shoe; 
however,  the  stress  in  these  areas  is  significantly  less  than  the  fracture  strength. 
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Figure  36  -  Maximum  Principal  Stress  for  394  m/s,  Coating,  Initial  Ambient 

Temperature 

The  effect  of  running  the  simulation  with  the  temperature  profile  can  be  seen  in 
Figure  37  where  the  highest  temperatures  are  found  toward  the  middle  of  the  shoe  at  x  = 
9.84  cm. 
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Figure  37  -  Temperature  and  Maximum  Principal  Stress  for  394  m/s.  Coating, 


Initial  Temperature  Profile 

Recall  that  the  temperature  profile  input  into  CTH  (Figure  32)  had  a  maximum 
temperature  of  525  K.  As  one  can  see  from  Figure  37,  the  effect  of  plasticity  raised  the 
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shoe  temperature  approximately  10  degrees.  The  maximum  prineipal  stress  deereased, 
yet  it  beeame  more  uniform  through  the  thickness. 


Flat  Rail,  798  m/s,  Coating 


As  expected,  increasing  the  horizontal  velocity,  and  therefore  the  vertical  velocity 
also  increased  the  resultant  temperature  in  addition  to  moving  the  area  of  high 
temperature  towards  the  front  of  the  shoe.  Figure  38  and  Figure  39  shows  a  higher 
temperature  initially  (at  1  microsecond  versus  10  microseconds)  of  approximately  422  K 
and  a  lower  maximum  principal  stress,  yet  these  values  are  both  still  well  below  what  is 
necessary  to  characterize  wear. 
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Figure  38  -  Temperature  for  798  m/s,  Coating,  Initial  Ambient  Temperature 
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Figure  39  -  Maximum  Principal  Stress  for  798  m/s,  Coating,  Initial  Ambient 

Temperature 

The  effect  of  the  temperature  profile  on  temperature  and  maximum  principal 
stress  can  be  seen  in  Figure  40.  Once  again  the  effect  of  the  initial  temperature  profile 
(Figure  33)  slightly  increased  the  temperature  and  thickened  the  portion  of  the  shoe 
experiencing  the  highest  value  of  maximum  principal  stress. 
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Figure  40  -  Temperature  and  Maximum  Principal  Stress  for  798  m/s.  Coating, 


Initial  Temperature  Profile 
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Flat  Rail,  394  m/s,  No  Coating 


The  effects  of  running  a  sled  directly  on  the  rail,  in  the  absence  of  a  coating,  will 
now  be  discussed.  Figure  41  shows  the  resulting  temperature  field  and  one-dimensional 
plot  of  temperature  through  the  thickness  of  the  shoe  at  the  end  of  10  microseconds  at  the 
X  position  of  20.5  cm. 
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Figure  41  -  Temperature  for  394  m/s,  No  Coating,  Initial  Ambient  Temperature 

The  effect  of  the  initial  temperature  profile  is  nearly  negligible  as  depicted  in 
Figure  42,  with  the  slight  distinction  between  the  far  field  temperature  of  550  K  versus 
300  K 
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Figure  42  -  Temperature  for  394  m/s,  No  Coating,  Initial  Temperature  Profile 


Flat  Rail,  798  m/s,  No  Coating 


Increasing  the  velocity  to  798  m/s  resulted  in  temperatures  that  can  be  seen  in 
Figure  43.  These  are  nearly  identical  to  the  394  m/s  case. 
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Figure  43  -  Temperature  for  798  m/s.  No  Coating,  Initial  Ambient  Temperature 

Again  the  effects  of  the  initial  temperature  profile  resulted  in  nearly  identical  results  for 
the  394  m/s  case  with  an  initial  temperature  profile  in  that  the  maximum  temperature 
increased  only  slightly. 
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Figure  44  -  Temperature  for  798  m/s,  No  Coating,  Initial  Temperature  Profile 


Asperity,  394  m/s,  Coating 

Cinnamon  [7:9-32]  used  an  asperity  to  generate  gouging  in  the  rail.  Asperities 
were  discussed  earlier  where  it  was  stated  that  they  exist  because  of  the  rail  seams  and  the 
inherent  manufacturing  tolerances.  His  study  conducted  on  asperity  geometry  (acute  face 
angle  and  height)  was  conducted  at  a  sled  velocity  of  3000  m/s  and  a  vertical  velocity  of 
-1  m/s.  He  concluded  that  the  highest  asperity  that  would  not  generate  a  gouge,  for  this 
velocity,  was  0.03  cm  with  a  face  angle  of  10  degrees.  This  particular  geometry,  well 
within  the  current  HHSTT  rail  specifications,  was  used  in  this  research  to  initiate  wear. 
Figure  45  shows  the  asperity  geometry. 
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Figure  45  -  Asperity  Geometry 

In  the  first  case  with  an  asperity,  the  speed  of  394  m/s  was  used  with  the  epoxy 
coating.  Figure  46  shows  a  high  temperature  spot  near  the  x  position  of  20.75  cm.  The 
graph  on  the  right  indicates  that  the  highest  temperature  at  this  position  is  very  nearly 
1500  K. 
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Temperature 
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Figure  46  -  Temperature  for  394  m/s,  Asperity,  Coating,  Initial  Ambient 

Temperature 

Figure  47  displays  that  the  maximum  principal  stress  at  an  x  position  of  20.72  cm 
is  significantly  greater  than  the  fracture  strength  of  VascoMax  300.  This,  according  to 
the  earlier  discussion,  is  an  indication  of  wear.  The  amount  of  wear  could  be  considered 
by  taking  the  thickness  of  the  shoe  that  is  above  the  fracture  strength  and  subtracting  off 
the  height  of  the  asperity  and  the  coating  (0.04  cm);  in  this  case  it  would  be 
approximately  0.08  cm. 
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Figure  47  -  Maximum  Principal  Stress  for  394  m/s,  Asperity,  Coating,  Initial 

Ambient  Temperature 

The  effects  of  applying  the  initial  temperature  profile  due  to  frictional  heating  in  this  case 
appear  to  be  negligible.  A  quick  comparison  of  Figure  46  with  Figure  48  and  Figure  47 
with  Figure  49  indicates  that  the  maximum  temperatures  and  maximum  principal  stresses 
are  nearly  equivalent. 
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Figure  48  -  Temperature  for  394  m/s.  Asperity,  Coating,  Initial  Temperature  Profile 
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Figure  49  -  Maximum  Principal  for  394  m/s,  Asperity,  Coating,  Initial  Temperature 


Profile 


Asperity,  798  m/s,  Coating 


This  case  is  identical  to  the  previous  case  with  the  exception  of  an  increase  in 
vertical  and  horizontal  velocities.  The  vertical  velocity  was  increased  to  -1.98  m/s  and 
the  horizontal  velocity  to  798  m/s  in  accordance  with  the  DADS  point  selected  earlier  in 
Figure  31.  As  Figure  50  indicates  the  temperature  nearly  reaches  1580  K. 
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Figure  50  -  Temperature  for  798  m/s.  Asperity,  Coating,  Initial  Ambient 


Temperature 


74 


Figure  5 1  shows  the  maximum  principal  stress  for  this  case,  and  like  the  slower 
velocity,  wear  is  also  evident  as  the  stress  exceeds  the  fracture  strength  of  the  material. 
However  this  case  reaches  an  even  higher  value  than  the  slower  velocity,  6.7  GPa  versus 
6.0  GPa.  As  a  result  the  wear  also  increases  by  a  thickness  of  0.05  cm  to  0.13  cm. 
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Figure  51  -  Maximum  Principal  Stress  for  798  m/s,  Asperity,  Coating,  Initial 

Ambient  Temperature 

The  effect  of  using  the  frictional  heat  temperature  profile  is  displayed  in  Figure  52  and 
Figure  53.  The  temperature  and  maximum  value  of  the  maximum  principal  stress  were 
significantly  reduced,  however  the  wear  thickness  increased  to  0.14  cm. 
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Figure  52  -  Temperature  for  798  m/s,  Asperity,  Coating,  Initial  Temperature  Profile 
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Figure  53  -  Maximum  Principal  Stress  for  798  m/s.  Asperity,  Coating,  Initial 


Temperature  Profile 


Asperity,  394  m/s,  No  Coating 

This  next  case  involves  the  slower  speed  of  394  m/s  case  with  an  asperity, 
however  there  is  no  epoxy  coating.  The  results  are  shown  in  Figure  54  and  Figure  55. 
Figure  54  indicates  that  the  maximum  temperature  is  approximately  1520  K,  while  Figure 
55  shows  a  wear  amount  of  0.16  cm  thick  resulting  from  the  maximum  principal  stress. 
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Just  as  in  the  case  with  this  velocity  with  a  coating,  there  is  no  appreciable  difference  at 
the  boundary  resulting  from  the  initial  temperature  profile  due  to  frictional  heating.  It 
can  be  seen  in  comparing  Figure  54  with  Figure  56  that  the  temperature  away  from  the 
boundary  is  elevated  due  to  the  profde,  but  this  appears  to  not  affect  the  wear  or 
maximum  temperature. 
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Figure  54  -  Temperature  for  394  m/s,  Asperity,  No  Coating,  Initial  Ambient 

Temperature 
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Figure  55  -  Maximum  Principal  Stress  for  394  m/s.  Asperity,  No  Coating,  Initial 


Ambient  Temperature 
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Figure  56  and  Figure  57  are  the  results  for  temperature  and  maximum  prineipal  stress 
when  an  initial  temperature  profile  is  used. 
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Figure  56  -  Temperature  for  394  m/s,  Asperity,  No  Coating,  Initial  Temperature 


Profile 


5  .  0 

4 . 5 
4 . 0 

3 . 5 
3 . 0 


>  1.5 


1  .  0 
0 . 5 
0  .  0 
-0  .  5 
-1  .0 


Maximum  Principal  Stress 


20.4  20.5  20.6  20.7  20.8  20.9  21.0 

I  2DR  X  (cm) 

-0.771  m/s,  394  m/sec,  no  coating,  asperity,  .002  cm  mesh,  TO=Prcfile 


M. 


BLLEVS  2/12/07  18:47:09  CTH 


6831  Time=1.00005x10"' 


(£1: 


I  2DR  Block  1  Y  Position  (10" ’em)  X=2.0725xt0’ 

-0.771  m/s,  394  m/sec,  no  coating,  asperity,  .002  cm  mesh,  T0=Profile 
BLLEVS  2/12/07  18:47:09  CTH  6831  Tlme=1.00003x10"5 


Figure  57  -  Maximum  Principal  Stress  for  394  m/s.  Asperity,  No  Coating,  Initial 


Temperature  Profile 
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Asperity,  798  m/s,  No  Coating 


The  final  case  in  this  combination  of  scenarios  is  the  798  m/s  case  with  the  shoe 
running  into  an  asperity  to  initiate  wear  and  no  coating  on  the  rail.  Figure  58  shows  that 
the  temperature  reaches  a  value  of  approximately  1590  K,  the  highest  temperature  of  all 
the  scenarios.  And  the  wear  as  indicated  by  the  maximum  principal  stress  being  greater 
than  the  fracture  stress  shows  a  value  of  0.23  cm  (Figure  59). 
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Figure  58  -  Temperature  for  798  m/s,  Asperity,  No  Coating,  Initial  Ambient 

Temperature 
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Figure  59  -  Maximum  Principal  Stress  for  798  m/s,  Asperity,  No  Coating,  Initial 

Ambient  Temperature 

The  temperature  profile  from  frietional  heating  affeeted  the  maximum  temperature  by 
redueing  it  by  approximately  100  K  to  1500  K,  yet  it  inereased  the  maximum  prineipal 
stress  to  6.4  GPa  and  inereased  the  wear  to  0.24  em.  These  effects  can  be  seen  in  Figure 
60  and  Figure  61. 
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Figure  60  -  Temperature  for  798  m/s.  Asperity,  No  Coating,  Initial  Temperature 
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Figure  61  -  Maximum  Principal  Stress  for  798  m/s,  Asperity,  No  Coating,  Initial 

Temperature  Profile 


Asperity,  1,411  m/s.  No  Coating 

One  final  case  was  conducted  to  determine  if  melt  would  occur  if  the  horizontal 
velocity  were  increased.  This  case  was  simulated  at  ambient  temperature  only.  Figure  62 
and  Figure  63  respectively  show  the  results  for  temperature  and  maximum  principal 
stress.  The  melt  temperature  was  not  reached;  in  fact  the  maximum  temperature  was 
1540  K,  which  was  less  than  the  similar  case  run  at  798  m/s.  On  the  other  hand,  the 
depth  into  the  shoe  that  the  maximum  principal  stress  exceeded  the  fracture  strength  was 
significantly  higher  than  all  of  the  preceding  cases.  That  thickness  was  0.37  cm.  Figure 
62  also  shows  behavior  similar  to  gouging  in  that  a  mass  of  material  is  flowing  toward 
the  rear  of  the  shoe,  however  there  is  no  mixing  with  the  rail  that  would  indicate  gouging. 
If  the  simulation  were  to  continue,  it  would  be  expected  that  this  mass  would  generate  a 
wear  particle. 
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Figure  62  -  Temperature  for  1,411  m/s,  Asperity,  No  Coating,  Initial  Ambient 
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Figure  63  -  Maximum  Principal  Stress  for  1,411  m/s.  Asperity,  No  Coating,  Initial 

Ambient  Temperature 

Table  7  is  a  list  of  the  results  from  the  17  numerical  cases  simulated  in  CTH. 
Listed  are  the  maximum  temperatures,  maximum  value  for  the  maximum  principal 
stresses  and  the  corresponding  wear  amounts,  if  any  occurred.  It  should  be  noted  that  the 
wear  values  listed  in  Table  7  are  not  material  removal  amounts.  They  simply 
characterized  the  depth  at  which  the  maximum  principle  stress  exceeded  the  fracture 
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strength.  Therefore,  these  values  cannot  be  compared  to  the  values  in  Table  6,  which 
represented  a  uniform  amount  of  material  removed  from  the  bottom  of  the  shoe. 

Table  7  -  Numerical  Wear  Results  from  CTH 


Max 

Max  Melt  Principal 


Temperature  Wear  Stress  Wear 

(K) _ (cm)  (GPa)  (cm) 


Coating 

Ambient 

329 

- 

1.4 

- 

394  m/s 

Profile 

533 

- 

0.9 

- 

No 

Ambient 

1500 

- 

- 

- 

Flat 

Coating 

Profile 

1510 

- 

- 

- 

Rail 

Coating 

Ambient 

425 

- 

1.25 

- 

798  m/s 

Profile 

775 

- 

1.25 

- 

No 

Ambient 

1520 

- 

- 

- 

Coating 

Profile 

1500 

- 

- 

- 

Coating 

Ambient 

1500 

- 

6.0 

0.08 

394  m/s 

Profile 

1500 

- 

6.0 

0.08 

No 

Ambient 

1520 

- 

6.2 

0.16 

Asperity 

Coating 

Profile 

1520 

- 

6.2 

0.16 

Coating 

Ambient 

1580 

- 

6.7 

0.13 

798  m/s 

Profile 

920 

- 

5.2 

0.14 

No 

Ambient 

1590 

- 

6.2 

0.23 

Coating 

Profile 

1500 

- 

6.4 

0.24 

Asperity 

1411  m/s 

No 

Coating 

Ambient 

1500 

- 

5.0 

0.37 
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IV.  Conclusion 


The  first  and  most  obvious  conclusion  is  that  melt  wear,  for  the  1 7  cases  studied, 
was  never  generated  in  CTH  because  the  shoe  material,  VascoMax  300,  never  reached  its 
melt  temperature  of  1685  K.  It  was  originally  thought  that  the  application  of  a  frictional 
heating  would  aid  the  shoe  in  reaching  elevated  temperatures  when  in  fact  it  did  very 
little  at  the  interface.  This  is  most  likely  due  to  the  fact  that  the  initial  temperatures 
changed  the  shoe  properties,  i.e.  thermal  softening,  enough  so  that  less  heat  was 
generated  due  to  plasticity.  In  other  words,  because  the  material  was  thermally  softened 
at  the  beginning  of  the  simulation,  less  energy  was  required  to  plastically  deform  the 
material,  and  therefore  less  energy  in  the  form  of  heat  was  generated. 

The  second  and  also  obvious  conclusion  is  that  the  fiat  rail  simulations  neither 
generated  melt  wear  nor  did  it  generate  wear  in  terms  of  material  failure.  In  fact  it  might 
be  fair  to  suggest  that  the  values  for  maximum  principal  stress  acquired  from  all  of  the 
simulations  on  the  fiat  rail  were  of  little  use.  In  terms  of  CTH,  sliding  the  shoe  on  the  flat 
rail  results  in  producing  mostly  shearing  stresses.  Looking  back  at  Figure  36  and  Figure 
39  shows  that  the  maximum  principal  stress  values  are  not  even  at  the  interface.  It  is  not 
until  there  is  some  type  of  impact,  such  as  with  an  asperity,  that  reasonable  values  for 
principal  stresses  contribute  to  wear. 

In  reality,  using  an  asperity  in  CTH  is  a  matter  of  practicality.  At  the  track  there 
are  rail  seams  approximately  every  15  m  (50  ft).  One  could  statistically  determine  how 
many  asperities  a  sled  would  encounter  at  a  given  speed  if  only  25%  of  these  seams,  for 
example,  had  a  wear  producing  asperity  (still  yet  to  be  determined).  A  technique  that 
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estimates  the  number  of  these  asperities  eneountered  for  a  given  run  eould  be  employed 
to  determine  the  total  wear  for  that  run. 

Previous  researeh  determined  that  gouging  is  initiated  by  asperities.  Keeping  the 
asperities  below  the  eritieal  values  for  gouging,  as  determined  by  Cinnamon,  may  not  be 
enough  to  eliminate  wear  generation  at  those  asperities  as  mentioned  in  the  introduetion. 

The  presenee  of  the  epoxy  eoating,  espeeially  in  the  eases  without  the  asperity, 
signifieantly  reduees  the  interfaee  temperature.  This  oeeurs  beeause  the  yield  strength  of 
the  epoxy  is  lower  than  the  yield  strength  of  either  the  steel  rail  or  the  VaseoMax  300 
shoe.  The  eells  at  the  interfaee  are  therefore  weaker  with  the  epoxy  and  require  less 
energy  to  meet  the  no-slide  eriteria  mentioned  in  ehapter  two.  In  effeet,  adding  the  epoxy 
lowered  the  eoeffieient  of  frietion.  In  the  ease  with  the  asperity,  the  eoating  also  reduees 
the  interfaee  temperature  but  to  a  mueh  lesser  degree  than  the  flat  rail. 

Finally,  the  values  for  wear  that  were  aehieved  using  the  asperity,  0.03  em  in 
height  and  10  degrees,  are  not  eonsistent  with  wear  values  produeed  by  the  DADS 
analytieal  integration  teehnique,  nor  are  they  realistie  for  a  1 0  mieroseeond  period  of 
time.  However  it  is  eoneeivable,  with  the  use  of  the  asperity,  perhaps  even  a  smaller  one, 
that  these  simulations  eould  eneompass  a  larger  duration  of  the  run  than  the  10 
mieroseeonds  in  real  time,  in  whieh  ease  the  wear  in  Table  7  would  be  more  reasonable. 

Additionally,  these  wear  values  were  measured  at  one  loeation  near  the  front  of 
the  shoe,  more  speeifieally,  at  the  eontaet  point  of  the  shoe  and  the  asperity.  This 
measurement  of  wear  depth  should  be  distributed  over  the  entire  length  of  the  shoe, 
whieh  would  reduee  the  thiekness  to  a  more  realistie  value. 
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As  demonstrated,  CTH  ean  be  used  to  show  wear  due  to  material  failure,  however 
quantifying  the  wear  in  relation  to  real-time  is  diffieult  beeause  the  simulations  run  for 
sueh  a  short  period  of  time.  The  author  is  eonfident  that  future  researeh  in  this  area  ean 
result  in  realistie  and  quantifiable  wear  values  due  to  material  failure  as  well  as  melting. 

Future  Research 

It  has  been  eonveyed  that  wear  is  a  eomplieated  problem  to  understand  and 
therefore  difficult  to  model.  Areas  of  future  research  should  include  attempting  to 
correlate  wear  predicted  by  the  analytical  model  and  wear  generated  by  CTH.  It  may  be 
reasonable  to  adjusting  the  parameters  studied  here  to  a  practical  dimension,  namely  the 
asperity  height,  to  find  the  values  which  produce  wear  in  amounts  that  match  the 
analytical  model  and  even  real  track  samples.  By  carefully  choosing  more  scenarios  in 
combinations  of  vertical  and  horizontal  velocities  (in  accordance  with  DADS),  asperity 
height  and  angle,  initial  temperature,  and  coating,  it  is  conceivable  that  a  function  could 
be  established  to  accurately  determine  wear  for  the  entire  sled  run. 

One  of  the  big  assumptions  made  in  this  research  was  with  respect  to  the  mass  of 
the  sled.  For  simplification,  it  was  assumed  that  the  entire  sled  mass  could  be  evenly 
divided  among  the  four  shoes,  and  in  addition,  that  mass  would  be  constant  through  the 
run.  In  fact,  the  mass  of  the  model  used  in  CTH  was  never  varied.  All  simulations, 
regardless  of  velocity,  used  one  quarter  of  the  entire  sled  mass.  A  quick  check  was 
performed,  after  the  fact,  to  see  if  the  change  in  momentum  over  time  would  be  equal  to 
the  vertical  force  supplied  by  DADS  using  the  equation 
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(28) 


^  t  —t 
‘2  ‘1 

The  results  from  this  proeedure  were  revealing.  Orders  of  magnitude  separate  the  two 
normal  foree  values  in  some  intervals,  while  others  eorrelated  very  well.  In  reality  the 
sled  is  pitehing  and  rolling,  and  this  ealculation  eonfirms  it.  These  effeets  ean  aetually  be 
seen  in  the  DADS  file  if  all  four  shoes  are  examined  together.  As  the  forees  inerease  on 
the  front  shoes,  the  forces  on  the  rear  shoes  decrease  indicating  the  sled  might  be  pitching 
forward.  This  suggests  that  the  mass  as  well  as  the  vertical  and  horizontal  velocities  need 
to  be  unique  input  parameters  for  each  CTH  simulation. 

Another  area  for  future  research  concentration,  as  alluded  to  in  the  conclusions,  is 
attempting  to  statistically  characterize  the  entire  sled  run  in  terms  of  shoe/rail  contacts. 
Figure  15  shows  these  contacts  over  a  brief  time  period.  If  these  events  could  be 
examined  for  the  entire  run  and  categorized  by  force,  duration/distance  traveled  and/or 
velocity  and  quantified  in  a  histogram  of  sorts,  then  selected  CTH  simulations  could  be 
executed  that  adequately  represent  each  category  of  contact.  The  results  from  that  run 
could  then  be  multiplied  by  the  number  of  contacts  and  summed  for  the  entire  run. 

DADS  does  not  currently  perform  such  calculations,  however  an  algorithm  could  be 
easily  generated  in  MATLAB. 

With  regard  to  the  analytical  model  and  the  development  of  kA  from  equation 
(17),  a  refined  range  of  kA  as  a  function  of  Pv  should  be  examined.  As  mentioned,  the 
range  of  Pv  for  the  shoe  was  0  to  35  GPa  m/s  however  almost  96%  of  the  Pv  values  for 
the  shoe  are  less  than  4.2  GPa  m/s.  The  average  Pv  value  for  the  entire  run  was  0.57  GPa 
m/s  which  is  not  as  telling  considering  all  of  the  zero  values  (when  there  is  not  contact 


87 


between  the  shoe  and  the  rail)  ineluded  in  this  average.  The  emphasis  however  is  on 
refining  the  range  on  whieh  the  kA  values  were  established.  This  would  give  a  better 
approximation  of  the  analytieal  wear.  On  the  same  note,  the  funetion  for  the  eoeffieient 
of  frietion  may  need  to  be  revisited  to  reduee  its  value  in  the  lower  Pv  range  (see  Figure 
20). 

The  models  used  for  the  DADS  integration  teehnique  were  based  on  the  pin-on- 
disk  geometry.  These  models  may  need  to  be  modified  and  adapted  to  the  shoe  and  rail 
system.  More  specifically,  the  pin  radius,  ro,  doesn’t  easily  translate  to  the  shoe.  This 
research  used  a  value  of  5.08  cm  (2  in),  half  the  width  of  the  shoe.  A  more  appropriate 
value,  based  on  the  ratio  of  areas,  would  have  been  8.10  cm  (2.19  in).  Also,  as 
previously  mentioned,  a  two  dimensional  heat  transfer  analysis  may  prove  to  be 
beneficial  as  well  as  considering  the  effects  of  using  material  properties  that  vary  with 
changing  temperature. 

With  regard  to  melt  wear,  additional  work  should  be  explored  concerning  the 
concept  of  flash  temperatures  at  the  interface  of  the  shoe  and  the  rail.  This  is  the  extreme 
temperature  that  occurs  on  the  microstructure  level  where  the  surface  asperities  contact. 
The  melt  equation  along  with  the  a  variable  for  heat  distribution  between  the  shoe  and 
rail  should  be  more  closely  examined  and  better  understood. 

For  the  cases  studied,  CTH  did  not  produce  temperatures  above  the  melting  point 
for  VascoMax  300,  suggesting  that  melt  wear  did  not  occur  for  the  cases  studied. 
Additionally,  the  heat  transfer  model  did  not  achieve  melting  either,  at  least  up  to  1,41 1 
m/s.  On  the  other  hand,  the  Lim  and  Ashby  wear  model  produces  melt  wear  for  the 
entire  sled  run  (see  Figure  21)  based  on  equation  (11).  This  could  be  that  both  CTH  and 
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the  heat  transfer  models’  eell  sizes  and  finite  differencing  meshes  are  not  small  enough  to 
detect  these  high  flash  temperatures  that  occur  very  close  to  the  boundary. 

Oxidation-dominated  wear  was  briefly  addressed  in  this  paper.  Recall  that  the 
values  were  unrealistically  large  when  computed  for  the  Montgomery  data.  This 
mechanism  of  wear  may  be  a  significant  contributor  to  the  overall  amount  of  wear  that 
takes  place.  In  chapter  two,  an  equation  representing  a  model  for  oxidation  wear  was 
given  for  the  pin-on-disk  scenario,  however  it  was  neither  explained  nor  examined 
beyond  calculating  oxidation  wear  for  the  DADS  data. 

Finally,  wear,  as  mentioned  several  times,  is  complex.  The  scenario  at  Holloman 
adds  to  the  complexity  due  to  the  speed  and  bouncing  effects.  This  will  not  be  a  problem 
that  is  easily  modeled  in  any  fashion.  Reliance  on  empirical  techniques  will,  without  a 
doubt,  prove  to  be  advantageous  in  this  effort.  In  light  of  that,  it  may  be  beneficial  to 
construct  a  test  setup  similar  to  the  pin-on-disk  and  connect  the  pin  to  an  actuator  that  can 
execute  an  input  function  that  mimics  a  particular  DADS  file  in  terms  of  vertical 
velocities  and  forces. 


89 


Appendix  A, 


MATLAB  Code 


9-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-9-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-9-9-5-9- 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Script  Wear2.m 

%  Capt  Greg  Cameron 
%  Thesis  -  Winter  2007 

%  This  script  reads  in  the  DADS  file  containing  time(t),  velocity (v), 

%  force (f)  and  vertical  velocity (vv) .  It  also  asks  the  user  for  a 
%  velocity  to  calculate  frictional  heating  up  to  that  input  value. 

%  Then  the  script: 

%  -  Calculates  distance  and  acceleration  and  plots 

%  along  with  velocity,  vertical  velocity,  force  (Thesis  Figure  14-15) 
%  -  Plots  points  vertical  and  horizontal  velocities 

%  vs  time  used  in  CTH  runs  (Thesis  Figure  31) 

%  -  Calculates  coefficient  of  friction  (Thesis  Figure  20) 

%  -  Calculates  wear  from  Dim  and  Ashby  equations  for 
%  melt  and  "mechanical"  wear  from  kA  found  in 

%  montgomery.m  script  using  Archard's  equation  (Thesis  Figure  21) 

%  -  Calls  Melt2.m  to  plot  melt  wear  vs  Dim  and 

%  Ashby  model  (Thesis  Figure  27-28) 

%  -  Calculates  the  melt  and  mechanical  wear  thickness 
%  by  integrating  as  well  as  the  heat  flux  (Thesis  Figure  22) 

%  -  Calls  Heat2.m  to  solve  differential  heat  equation 
%  to  find  temperature  profile  in  the  shoe  from  the 

%  heat  flux  boundary  condition  (Thesis  Figure  32-33) 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

g, 

o 

clc 

clear  all 
close  all 
warning  off  all 
format  long 

%  Prompt  user  to  input  ending  velocity  in  m/ s 
velocity  =  input (' Input  Velocity  in  m/s:\n\n'); 


%  Definition  of  Constants 
Width  =4;  % 
Length  =8;  % 
An  =  Width*Length;  % 
rO  =  2;  % 
a  =  .006993;  %. 014105;  % 
H  =  1.45e5;  % 
beta  =6;  % 
TO  =  535;  %522  % 
Tstar  =  264;  % 
Tm  =  3033;  % 
L  =  3.05e5;  % 
lb  =  .04;  % 
KA  mild  =  5e-5;  % 
KA_severe  =  5e-3;  % 
K  =  0.5;  % 
mass  DADS  =1854/4;  % 


Shoe  Width 

Shoe  Length 

Nominal  Contact  Area 

Radius  of  Nominal  Contact  Area 

Thermal  Diffusivity  (Vascomax) 

VascoMax  Hardness 

Numerical  Constant 

Sink  Temperature  (522  Lim/Ashby) 

=aH/K 

Melting  Point  for  VascoMax  (1685  K) 
Latent  heat  of  melting  for  steel 

Archard's  Coefficient  for  mild  wear 
Archard's  Coefficient  for  severe  wear 
Fraction  of  heat  transfer 
From  Holloman  on  18  Jan  07 


in 

in 

in''2 


in 

in'^2  /  s 
psi 


R 

R 

R 

psi 

in 


Ibm 
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%f ile= [ ' I : \My  Documents\Thesis\DADS . txt ' ] ; 

file=  [  '  z  :  \Greg\Work;\Air  Force \AFIT\Thesis\MATLAB\ DADS 3x  .  txt '  ]  ; 

%  time (s) , velocity (in/s) , force (Ibf) , vertical  velocity  (in/s) 
[t,v,f,vv]  =  textread (file, ' %f  %f  %f  %f'); 

%  Convert  v  and  vv  to  ft/s 
v=v/ 12 ; 
vv=vv/ 12 ; 
n  =  size  (t, 1 )  ; 


%Calculate  the  distance  for  each  step  by  integrating  the  velocity 

d(l)=0; 

accel ( 1 ) =0 ; 

for  i=2 : n 

d  (i)  =d  (i-1)  +  (v  (i) +v  (i-1)  )  *  (t  (i) -t  (i-1)  ) /2;  %In  feet 

accel (i)=(v(i)-v(i-l) ) / (t(i)-t(i-l) ) ; 

end 


%  Convert  to  SI  units 
d_SI=d*. 30498; 
v_SI=v*. 30498; 
vv  SI=vv* . 304 98 ; 
accel_SI=accel* .30498; 
f  SI=f*4.448; 


ooooooooooooooooooooooooooooooooooooooooo 

%  Plot  dynamics  of  the  entire  sled  run  % 

ooooooooooooooooooooooooooooooooooooooooo 


figure  ( 1 ) 

set (gcf , ' Color ' , 

subplot  (4,2,1); 


subplot  (4,2,2)  ; 


subplot ( 4 , 2 , 3 : 4 ) 


subplot ( 4 , 2 , 5 : 6 ) 


subplot  ( 4 , 2 , 7  :  8 ) 


%Thesis  figure  14-15 

'w'  ) 

plot (t, d_SI , ' b ' ) ; ylabel ( ' Distance  (m) ' ) ; 

xlabel ( ' Time  ( s ) ' ) ; 
xlim([0  2.5] ) ;%ylim( [4767  4844]); 
plot ( t , v_SI ,  ' b ' ) ; ylabel (' Velocity  (m/s )  ' )  ; 

xlabel ( ' Time  ( s ) ' ) ; 

xlim([0  2 . 5 ] ) ; ylim ( [ 0  3000]); 

%ylim( [6292  6296] ) ; 

;  plot  ( t ,  accel_SI ,  '  b  '  )  ;  ylabel  (  '  Acceleration  (m/s''2 )  '  )  ; 

xlabel ( ' Time  ( s ) ' ) ; 

xlim([0  2.5] ) ;ylim( [0  2438]); 

; plot (t, vv_SI ,' b '); ylabel ( 'Vertical  Velocity  (m/s)'); 

xlabel ( ' Time  ( s ) ' ) ; 
xlim([0  2.5] ) ;ylim( [-3  3] ) ; 
%xlim([.3  .32] ) ;ylim( [-1.5  3]); 

; plot (t, f_SI ,' b '); ylabel (' Vertical  Force  (N) ' ) ; 

xlabel ( ' Time  ( s ) ' ) ; 
xlim ( [ 0  2.5]); 

%xlim([.3  .32]); 


5-9-5-9-5-9-5-Q-5-9-5-9-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-9-9-5-Q-5-9-5-9-5-9-9-Q-5-9-5-9-5-9-5-Q-5-9-5-9-5-9-5-Q-5-9-5-9-5- 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Plot  Velocity  and  vertical  velocity  vs  time  showing  points 
%  considered  in  CTH  Modeling 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

figure  (2)  %Thesis  Figure  31 

set (gcf ,  'Color',  'w') 
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o\o  o\o 


subplot (2,1,1); plot ( t , v_SI , 
hold  on 


' b ' ) ; ylabel ( ' Down 
xlabel ( ' Time 


plot (. 31 , 394 . 0634 ,' ro '); hold  on 
plot(. 6298, 797. 5336, 'ro');hold  on 
plotd. 1226,  1410.8598,  '  ro  '  ) 
text (.29,800, 'A') 
text ( . 61, 1200, 'B' ) 
text  (1 . 1,  1816,  'C  ) 


Track  Velocity  (m/s) ' ) 
(s)  ' ) ;%xlim( [ . 08  . 2 ] )  ; 


subplot (2,l,2);plot(t, vv_SI , ' b ' ) ; ylabel (' Vertical  Velocity  (m/s ) ' ) ; 

xlabel ( ' Time  ( s ) ' ) ; %xlim ( [ . 08  . 2 ] ) ; 


hold  on 

plot (. 31, -. 7714, ' ro ') ;hold  on 
plot(. 6298, -1.98330, 'ro');hold  on 
plotd. 1226, -3. 47610,  '  ro  '  ) 
text (.32,-2, ' A ' ) 
text ( . 66, -3 . 1 ,  ' B ' ) 
text(1.15,-4.45, 'C') 


%  Calculate  the  pressure,  pv,  Ftil,  vtil  and  pvtil  for  each  step 

p=f./An;  %pressure  (psi) 

pv=p. *v*le-6;  % (psi*ft/s*le-6) 

pv_SI=pv/le-6* . 30498*6 . 895e-6;  %GPa  m/s 

Ftil=f . /An. /H; 

vtil=v .*r0./a.*12; 

pvtil=Ftil . *vtil ; 

%  Calculate  coefficient  of  friction  and  plot  versus  pv 
mu=1.5704.* (pvtil) .^-.2299; 

C=min (mu) -0.02; 
mu=mu-C; 

Bravo=isnan (mu) ; 
for  z=l : n 

if  Bravo (z) ==1 
mu ( z ) =  0 ; 
end 

end 


figure (3)  %Thesis  Figure  20 

plot (pvtil, mu, 'b. ' ) 
ylim( [0  0.8] ) 

set (get (gca, ' XLabel ' ) , ' Interpreter ' , ' Latex ' ) 
xlabel ( '$\tilde{P}\tilde{v}$' ) 

ylabel (' Coefficient  of  Friction  (  \it\mu\rm) ') 
set (gcf , 'Color', 'w') 

%alpha=l .  /  (2+beta*  (pi/8  .  *vtil)  .  '^  ( .  5)  )  ; 
alpha=l./ (2  +  lb* (pi/8/a/rO . *v*12)  (.5) ) ; 

%  Calculate  Melt  Wear 

Wmelt= ( (Tm-TO) /Tstar*H/L/beta)  . /vtil . * (alpha. *mu. *Ftil . *vtil* .  .  . 

(Tstar*beta/ (Tm-TO) ) -1) *An*12; 

Wmelt  SI=Wmelt*53731 . 6;  %in  (mm''3/m) 

Gamma=isnan (Wmelt)  ; 
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Phi=isnan (Wmelt  SI); 
for  z=l : n 

if  Gamma (z)==l 
Wme 1 1 ( z ) =  0 ; 

end 

if  Phi ( z ) ==1 

Wmelt_SI (z) =0; 

end 

end 

%  Calculate  Mild  and  Severe  Plastic  Wear 
Wplas  mild=KA  mild*Ftil*An*12; 

Wplas  mild  SI=Wplas  mild*53731 . 6;  %in  (mm''3/m) 

Wplas  severe=KA  severe*Ftil*An*12; 

Wplas  severe  SI=Wplas  severe*53731 . 6;  %in  (mm''3/m) 

%  Determine  kA  based  on  experimental  Montgomery  Data 
for  j  =1 :  n 

if  pv ( j ) <2 . 00 

kA ( j ) =3 . 534e-3 ; 
elseif  pv(j)>=2.00 
if  pv(j)<=15.00 

kA ( j ) =8 . 1 60e-4 ; 
else 

kA ( j ) =1 . 303e-3; 

end 

end 

end 


kA=kA' ; 

W=kA. *Ftil*An*12; 

W_SI=W*53721 . 6; 

figure (4)  %Thesis  figure  21 

set (gcf , 'Color', 'w') 

subplot (2, 1, 1) ;plot (d_SI , W_SI , ' b ' , d_SI , Wmelt_SI , ' r ' ) ; 

ylabel('Wear  Rate  (mm^3/m) ') ;xlabel (' Distance  (m) ' ) ; 
legend ( ' \itw  k_A\rm' , ' \itw_m_e_l_t\rm' ) ; 
ylim([0  537 . 6] ) ;xlim( [0  3871]) 
subplot (2, 1,2), -plot (d_SI , W_SI , ' b ' , d_SI , Wmelt_SI , ' r ' ) ; 

ylabel('Wear  Rate  (mm''3/m)  ');  xlabel  (' Distance  (m)  '  )  ; 
legend ( ' \itw_k_A\rm' , ' \itw_m_e_l_t\rm' ) ; 
ylim([0  268.8] ) ;xlim( [575  625]) 

figure (5)  %Deleted  this  figure  from  Thesis 

set  (gcf,  'Color',  'w') 

plot (v_SI,W_SI, 'b' , v_SI, Wmelt_SI, ' r ' ) ; 

ylabel('Wear  Rate  (mm''3/m)  '  )  ; 

xlabel ( 'Velocity  (m/s)  '  )  ; 

legend ( ' \itw_k_A\rm' , ' \itw_m_e_l_t\rm' ) 

ylim([0  100] ) ;xlim( [350  450]) 

Mel'b2  %figure(6)  -  Thesis  figure  29-30 

%  Calculating  the  Wear  Volume  for  the  distance  traveled  by 
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%  integrating  wear  rate  as  a  function  of  distance  using  the 
%  trapezoid  rule 

WearP ( 1 ) =0 ; 

WearM ( 1 ) =0 ; 

Power ( 1 ) =0 ; 
i=2; 

%  Calculate  the  melt  and  mechanical  wear  and  heat  flux  up  to  the 
%  velocity  specified  by  the  user 

while  v(i)  <  velocity/ . 30498  %Convert  input  "velocity"  to  ft/s 

A=d ( i - 1 ) ; 

B=d ( i ) ; 

h=(B-A);  %the  distance  between  to  data  points 

WearP  (i)  =WearP  (i-1)  +  (W  (i-1) +W  (i)  )  *h/2;  %In  in^'B 

WearM ( i ) =WearM ( i-1 )  +  (Wmelt ( i-1 ) fWmelt ( i ) ) *h/2  ; 
if  mu ( i ) ==Inf 
mu ( i ) =  0 ; 

end 

Power (i) =K*mu (i) *f (i) *v (i) ;  %In  lbf*ft/s 

time ( i ) =t ( i ) ;  %For  plotting  purposes 

i=i+l ; 

end 

Flux=Power/An*144;  %In  Ibf/s/ft 

figure (7)  %Thesis  figure  22 

set(gcf, 'Color', 'w') 
subplot  (2 , 1 , 1 ) 

Flux_SI=Flux*14 .596; 
plot (time, Flux_SI , 'b') 
xlabel ( ' Time  ( s ) ' ) 
ylabel  (  '  Flux  (W/m''2 )  '  ) 
xlim ( [ 0  . 3 ] ) 

tf=t(i-l);  %Time  were  sled  reaches  input  speed 

hold  on 

%  Calculate  the  temperature  distribution  in  the  shoe 

Heat2 (TO, Flux, a, tf, velocity)  %figure(8)  -  Thesis  figure  32-33 

fprintf ( ' \nTotal  Distance  Traveled:  %6.3f  ft\n',  d(i-l)) 
fprintf ( '  or:  %6.3f  m\n',  d (i-1) *. 30498) 

fprintf (' This  corresponds  to  a  velocity  of:  %6.3f  m/s\n',... 
v(i-l) *.30498) 

fprintf ( '  or:  %6.3f  f t/s\n ' , v ( i-1 ) ) 

fprintf ('The  vertical  velocity  for  this  speed  is:  %6.3f  m/s\n',... 
vv (i-1) *.30498) 

fprintf  (' \nTotal  Plastic  Wear  for  This  Distance  is:  %12.9f  in''!',... 
WearP (i-1) ) 

fprintf (' \nThe  Reduction  in  Shoe  Height  is:  %12.9f  in ', WearP ( i-1 ) /An) 
fprintf (' \n  or:  %12.9f  cm\n ', WearP ( i-1 )/.. . 

An  *2.54) 

fprintf (' \nTotal  Melt  Wear  for  This  Distance  is:  %12.9f  in^3',... 

WearM ( i-1 ) ) 

fprintf (' \nThe  Reduction  in  Shoe  Height  is:  %12.9f  in',  WearM ( i-1 ) /An) 


%Total  Plastic  Wear  at  step  1  equals  0 
%Total  Mechanical  Wear  at  step  1  equals  0 
%Power  as  a  function  of  time 


94 


fprintf('\n  or:  %12 

WearM ( i-1 ) /An *2 . 54  ) 

fprintf ( ' \nThis  Occurs  at  Time  =  %12.9f  s\n\n', 
return 


9f  cm\n ' , . . . 
time (i-l) ) 
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%  Script  Melt2.m 

%  Capt  Greg  Cameron 
%  Thesis  -  Winter  2007 

%  This  script  is  called  from  Wear2.m  and  does  the  following: 

%  -  Creates  a  model  for  melt  wear  based  on  the  Lim  and  Ashby  model 

%  plots  the  DADS  data  that  was  evaluated  for  melt  wear  in  Wear2.m 

%  and  creates  a  video  file  to  shoe  in  3D  (Thesis  Figures  29-30) 

g, 

o 

9-9-5-9-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-9-Q-5-9-5-Q-5-9-5-Q-5-9-5-9-5-9-9-Q-5-9-5-9- 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Generating  Data  to  plot  contour  lines  for  Lim  and  Ashby  Model 

9-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-9-9-5-9-5-9-5-9-9-9-5-9-5-9-5-9-5-9-5-9-Q-9- 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

rO  model=2; 
lb_model= . 04  ; 
beta_model=6; 

An  model=pi*r0  model^2; 

X=logspace ( -2 , 8 , 1 00 ) ;  %This  is  Vtilde 
Y=logspace ( -5 , 5 , 1 00 ) ;  %This  is  Ftilde 

[XI, YI]  =  meshgrid (X, Y) ; 

vXY=XI *a* 12 /rO  model; 
muXY=1.5704.* (XI.*YI) (-.2299) ; 

%  C=min (min (muXY) ) -0 . 02  ; 

%  muXY=muXY- . 02; 

alphaXY=l . / (2  +  lb  model* (pi/8/a/rO  model . *vXY* 12 )  . ^  ( . 5 ) ) ; 

WmeltXY= ( (Tm-  ~  ~ 

TO) /Tstar*H/L/beta  model) . /XI . * (alphaXY. *muXY. *YI . *XI* (Tstar*beta  model 
/ (Tm-TO) ) -1) *An  model* 12; 


V=[le-9  le-8  le-7  le-6  le-5  le-4  le-3] ; 

figure (6)  %Thesis  Figure  29-30 

contour3 (XI , YI , WmeltXY, V) 
set (gcf , 'Color', 'w') 

set (get (gca, ' XLabel ' ) , ' Interpreter ' , ' Latex ' ) 
set (get (gca, ' XLabel '), 'String', '$\tilde{v}$') 
set (get (gca, ' YLabel ' ) , ' Interpreter ' , ' Latex ' ) 
set (get (gca, ' YLabel '), 'String', '$\tilde{F}$') 
set (get (gca, ' ZLabel ' ) , ' Interpreter ' , ' Latex ' ) 

set (get (gca,  ' ZLabel '),  'String',  '$\tilde{w}$',  ' Font Size ',16,  ' Rotation ', 0 

) 

zlim ( [ le-9  le-3 ] ) 
hold  on 


plot 3 (vtil , Ftil , Wmelt , ' b . ' ) 


set (gca, 'ZScale', 'log') 
set (gca, 'XScale', 'log') 
set (gca, 'YScale', 'log') 
set (gca, 'XTick', [.01  .1  1 
100000000] ) 

set  (gca,  'YTick',  [.00000001 
10  100  1000]) 

set (gca,  ' XMinorCrid '  ,  'off 


10  100  1000  10000 
.0000001  .000001 
, 'YMinorCrid' , 'off 


100000  1000000  10000000 
.00001  .0001  .001  .01  .1  1 
' , 'ZMinorCrid' , 'off ) 


96 


set (gca,  ' XMinorTick '  ,  ' of f '  ,  ' YMinorTick '  ,  ' of f  ,  ' ZMinorTick '  ,  ' of f ' ) 
ylim([le-7  le3]) 

1  =  0; 

J=90; 

set (gca, 'View' , [I  J] ) 

%return 

Az=0; 

El=90; 

KAz=0;  %Gain 

KE1=1;  %Gain 

mov  =  avif ile ( ' DADS_Melt . avi '  ,  ' EPS ' , 2 0 ) ; 

M(i)  =  getframe (gcf ) ; 
mov  =  addf rame (mov, M ( i ) ) ; 
for  i=l:150 

Az=Az-KAz ; 

E1=E1-KE1; 
if  El<61 
KAz=l; 

KE1=0; 

end 

if  Az<-109 
KAz=0; 

KE1=0; 

end 

Look= [Az  El ] ; 

set (gca, ' View ' , Look) ; 

M(i)  =  getframe (gcf ) ; 
mov  =  addf rame (mov, M ( i ))  ; 

end 

mov  =  close (mov) ; 
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function  []=Heat2 (Tf , Flux , alpha , tf , velocity) 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Capt  Greg  Cameron 
%  Thesis  Winter  -  2007 

%  This  program  solves  the  one-dimensional  heat  equation  as  called  from 
%  Wear2.m  and  calls  CTH_Prof ile . m  to  generate  the  points  needed  for 
%  creating  the  temperature  profile  in  CTH 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Tf  =  Sink  Temperature  of  the  sled  in  Rankin  of  the  sled  (Ambient) 

%  Flux  =  Flux  at  the  Boundary  as  a  function  of  time 
%  alpha  =  thermal  diffusivity  in  in^2/s 

%  tf  =  time  required  to  reach  input  velocity  in  seconds 

Flux=Flux* . 0014594 ;  %Convert  from  Ibf/s/ft  to  W/cm^2 

k=0. 306582;  %Thermal  conductivity  of  Vmax  in  W/cm/K 

alpha  =  alpha*2 . 54 ^2 ;  %Converting  to  cm''2/s  -  for  CTH 

Tf= (Tf-459 . 67-32) *5/9+273 . 15;  %Converting  to  Kelvin  -  for  CTH 

t0=0;  %Time  =  0 

dt=.0001;  %Time  Step 

depth=2.5;  %Depth  of  thermal  profile  into  shoe  (cm) 

dy=. 015; 

%  Check  Stability  as  determined  by  ratio  r 
r=alpha*dt/dy^2  ; 
fprintfCr  =  %6.3f\n',r) 
if  r>.5 

fprintf('r  should  be  less  than  0.5\n\n') 
return 

end 

%T(y,t)  =>  Temperature  as  a  function  of  position (y)  and  time(t) 
t=t0 : dt : tf ; 

N=size (t, 2 )  ; 

y=0 : dy : depth; 

M=size (y, 2 ) ; 

T(1:M, l)=Tf;  %Temp  just  over  the  interface  at  time  t=0  (Ambient) 

T(M, l:N)=Tf;  %Temp  far  away  from  the  interface  for  all  times  (Ambient) 

Tzero=T (2,1) +2*dy*Flux (1) /k; 
for  j=2:N  %time(t) 

T (1,  j)=T (1,  j-1) +r* (T (2,  j-1) -2*T (1,  j-1) +Tzero)  ; 
for  i=2:M-l  %position(y) 

T(i,  j)=T(i,  j-l)+r* (T(i  +  1, j-l)-2*T(i, j-l)+T(i-l, j-1) )  ; 

end 

Tzero=T (2, j ) +2*Flux ( j ) *dy/k; 

end 

%  This  goes  through  and  plots  the  Temperature  profile  into  the  shoe  as 
%  a  function  of  time.  Note  that  the  last  time  step  may  not  produce  the 
%  hottest  temperature  at  the  boundary  due  to  the  changing  flux, 
figure ( 8 ) 
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set (gcf , 'Color', 'w') 
for  k=l:l:j 

Temp=T ( : , k) ; 
plot (y ' , Temp) 

end 

hold  on 

T0=max (max (T) ) ;  %This  is  the  maximum  over  all  times  -  see  above 
xlim ( [ 0  . 5 ] ) ; 

%ylim( [Tf  TO] ) ; 
ylim( [250  800] ) ; 

YT=250:50:800; 
set (gca, 'YTick',YT) 

xlabel (' Vertical  Distance  into  the  Shoe  (cm) ') 
ylabel ( ' Temperature  (K) ' ) 

figure  ( 7 ) 

subplot (2 , 1,2) 

plot (t, T  (1,  :  )  )  ; 

xlabel ( ' Time  ( s ) ' ) 

ylabel ( ' Temperature  (K) ' ) 

xlim([0  .33]) 

%  Prepare  data  and  call  function  to  create  the  output  temp  profile 
%  for  input  into  CTH 

%  Temp  ends  up  being  the  last  column  in  T(y,t) 

CTH  Prof ile (y ', Temp)  %figure(8)  -Thesis  figure  32-33 

fprintf ( ' \nThe  hottest  temperature  reached  was:  %6.2f  K\n',  TO) 


99 


function  [ ] =CTH_Prof ile (y , T) 

9-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-9-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9- 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Capt  Greg  Cameron 
%  Thesis  -  Winter  2007 

%  This  file  will  analyze  the  last  temperature  vs  y  profile  from  Heat 

%  and  select  the  best  points  for  creating  a  temperature  profile  in  CTH 

0, 

o 

%  Inputs:  y  is  the  y  position 

%  T  is  the  temperature  at  that  y  position 

%  The  temperature  exponentially  decays  with  y  position  so  we  only  need 
%  to  take  points  that  change  significantly 

%  first  look  at  the  total  temperature  difference  from  y=0  to  y=endpoint 

%  y=0:100; 

%  T=200:-2:0; 

K=11604.5;  %Constant  for  converting  from  Kelvin  to  eV  for  CTH 

N=size (y, 1 ) ; 

T0=T  (1) ; 

Tf=T (N) ; 

[C, I ] =max (T) ;  %Find  the  Maximum  temperature 

DeltaT=C-Tf ;  %it  may  not  be  at  the  boundary 

%  First  Row  of  the  Temperature  Profile  for  CTH 
CTH(1,  :)  =  [y(l)  T(l)  T(l)/K]; 

Q=2;  %Now  fill  in  the  remaining  rows  until  we  get  to  Tf 

%  Take  points  at  1-Pct%  of  previous  Temperature  point. 

Pct=0.93;  %0nly  need  to  change  this  variable 
%(.93  for  394.1  m/s  and  798  m/s) 

D=l-Pct; 

Point=Pct*C; 
for  i=2:N 

if  T(i)<Point 

CTH(Q,  :)  =  [y(i)  T(i)  T(i)/K]; 

Q=Q+l; 

Pct=Pct-D; 

Point=Pct*C; 
if  Point<Tf+0.1 

while  T(i)>Tf+0.1 
i=i+l ; 

end 

CTH(Q,  :)  =  [y(i)  T(i)  T(i)/K]; 
break 

end 

end 

end 

figure (8)  %Thesis  figures  32  and  33 

plot (CTH (:, 1) ,CTH(:,2) , 'k. ') 
save  CTH. out  CTH  -ASCII 
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%  Script  Montgomery . m 

%  Capt  Greg  Cameron 
%  Thesis  -  Winter  2007 

%  This  script  reads  in  the  data  from  the  Montgomery  pin-on-disk 
%  experimentation  table  10.  Specifically  it  reads  in  velocity, 

%  pressure,  mu,  distance,  wear.  Then  the  script: 

%  -  Calculates  wear  from  the  Lim  and  Ashby  equations  for  melt 
%  wear  and  both  mild  and  severe  plastic  wear.  (Thesis  Figure  16) 

%  -  Calculates  kA  for  three  discrete  ranges  of  Pv  (Thesis  Figure  17) 

%  -  Recalculates  "mechanical"  wear  from  new  kA  (Thesis  Figure  18) 

%  -  Plot  experimental  coefficient  of  friction  with 
%  curve  fit  and  adjusts  the  curve  to  approach  0.02 

%  for  the  Pv  range  (Thesis  Figure  19) 

q, 

o 

%  Values  are  provided  in  English  as  well  as  SI  units.  Variables  with 
%  SI  units  are  denoted  with  the  SI  suffix  in  the  name. 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


clc 

clear  all 
close  all 


warning  off  all 
format  long 


%  Definition  of 
An  =  .00502655; 
rO  =  .04; 
a  =  .014105; 

H  =  1.45e5; 
beta  =  6; 

TO  =  535; 


Constants 

%Nominal  Contact  Area 

%Radius  of  Nominal  Contact  Area 

%Thermal  Diffusivity  (Steel) 

%Steel  Hardness 

%Numerical  Constant 

%Sink  Temperature  (522  Lim/Ashby) 


Tstar  =  400; 

Tm  =  3240; 

L  =  3.05e5; 
lb  =  .04; 

KA  mild  =  5e-5; 

KA  severe  =  5e-3; 


%=aH/K 

%Melting  Point  for  Steel  (1800  K) 
%Latent  heat  of  melting  for  steel 
%Linear  diffusion  distance  into  pin 
%Archard's  Coefficient  for  mild  wear 
%Archard's  Coefficient  for  severe  wear 


in''2 

in 

in^2/s 

psi 

R 

R 

R 

psi 

in 


%f ile= [ ' I : \My  Documents\Thesis\DADS . txt ' ] ; 

f ile= [ ' z : \Greg\Work\Air  Force\AFIT\Thesis\MATLAB\Mont5 . txt ' ] ; 


%  velocity(ft/s),  pressure  (psi )  ,  ( - )  ,  distance  (ft)  ,wear(in''3/ft*le-7) 
[v, p, mu, d, Wexp] =textread ( f ile,  ' %f  %f  %f  %f  %f'); 

Wexp=Wexp* le-7 ; 

Wexp  SI=Wexp*53731 . 6;  %Convert  to  SI  units  (mm'^3/m) 

n=size (v, 1 ) ; 

pv=p. *v*le-6;  % (psi*ft/s*le-6) 

pv_SI=pv/le-6* . 30498*6 . 895e-6;  % (CPA  m/s) 

Ftil=p/H; 

vtil=v .*r0./a.*12; 

Pvtil=Ftil . *vtil ; 
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%  Calculate  alpha  for  melt  wear  equation 
%alpha=l .  /  (2+beta*  (pi/8  .  *vtil)  .  .  5)  )  ; 

alpha=l./ (2  +  lb* (pi/8/a/rO . *v*12)  (.5) ) ; 

%  Calculate  Melt  Wear 

Wmelt= ( (Tm-TO) /Tstar*H/L/beta) . /vtil . * (alpha. *mu. *Ftil . *vtil* . . . 

(Tstar*beta/ (Tm-TO) ) -1) *An*12; 

Wme It  SI =Wme lt*53731 . 6; 

%  Calculate  Mild  Plastic  Wear 
Wplas  mild=KA  mild*Ftil*An* 12 ; 

Wplas  mild  SI=Wplas  mild*53731 . 6; 

%  Calculate  Severe  Plastic  Wear 
Wplas  severe=KA  severe*Ftil*An*12; 

Wplas  severe  SI=Wplas  severe*53731 . 6; 

%  Plot  the  3  wear  equations  with  experimental  versus  Pv 

figure (1)  %  Thesis  Figure  16 

%plot (pv, Wexp, ' b . ' , pv, Wplas_mild, ' c . ' , pv, Wplas_severe, ' g . ' , pv, . . . 

%  Wmelt, ' r . ' ) 

plot (pv  SI ,  Wexp  SI ,  ' b .  ' , pv  SI ,  Wplas  mild  SI ,  ' c .  ' , pv  SI ,  .  . . 

Wplas  severe  SI , ' g . ' , pv  SI , Wmelt  SI , ' r . ' ) 
legend (' Experimental Mild  Plastic Severe  Plastic Melt ' ) 

%ylim ( [ 0  3e-5 ] ) 
ylim([0  1.6]) 

%xlabel  (  ' \itPv\rm  (psi  ft/s  x  10''-''6)  ') 
xlabel ( ' \itPv\rm  (GPa  m/s)') 

%ylabel  (  '  Wear  Rate  (in''3/ft)') 
ylabel('Wear  Rate  (mm^'S/m)  ') 
set (gcf , 'Color', 'w') 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Solve  for  kA  using  empirical  data 

9-9-5-Q-5-9-5-9-5-9-5-Q-5-9-9-9-5-9-5-Q-5-9-5-9-9-9-5-Q-5-9-5-Q-9-9-5-Q-5-9-9-9-5-9-5-Q-5-9-9-Q-9-9-5-9-5-9-9-9-5-9-5-Q-5-9-5-Q-9-9-5-Q-5- 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

kA=H/An* (Wexp. /p) /12; 

%  This  is  for  kA  adjust3e.  If  kA  is  out  the  limits  set  by  Lim  and 
%  Ashby,  manually  bring  it  back  to  the  limit 
for  i=l:n 

kA_ad j ( i ) =kA ( i ) ; 
if  kA ( i ) >5e-3 

kA_adj (i)=5e-3; 

end 

end 

%  This  is  for  kA  averaged.  Define  3  arbitrary  regions  of  Pv,  average 

%  over  the  three  regions. 

kA  avg  1=0; 

kA  avg  2=0; 

kA  avg  3=0; 

countl=0 ; 

count2=0 ; 

count3=0 ; 

for  j=l:n 
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if  pv ( j ) <2 . 00 

kA  avg  l=kA  avg  1  +  kA  adj  (j)  ; 
count l=count 1+1 ; 
elseif  pv(j)>=2.00 
if  pv(j)<=15.00 

kA  avg  2=kA  avg  2  +  kA  adj  (j)  ; 
count2=count2+l ; 

else 

kA  avg  3=kA  avg  3+kA  adj (j) ; 
count3=count3+l ; 

end 

end 

end 

%  These  are  the  3  discrete  values  for  kA  as  a  function  of  Pv. 
for  m=l : n 

if  pv (m) <2.00 

kA  avg (m) =kA  avg  1 / count 1; 
elseif  pv(m)>=2.00 
if  pv (m) <=15 . 00 

kA  avg (m) =kA  avg  2/count2; 

else 

kA  avg (m) =kA  avg  3/count3; 

end 

end 

end 

kA  of  Pv  l=kA  avg  1/countl; 
kA  of  Pv  2=kA  avg  2/count2; 
kA  of  Pv  3=kA  avg  3/count3; 

%  Plot  the  3  series  in  solving  for  kA  from  experimental  data 
figure (2)  %  Thesis  Figure  17 

%plot (pv, kA, ' b . ' , pv, kA_adj , ' ro ' , pv, kA_avg, ' g . ' ) 
plot (pv_SI , kA, ' b . ' , pv_SI , kA_adj , ' ro ' , pv_SI , kA_avg, ' g . ' ) 

%xlabel ( ' \itPv\rm  (psi  ft/s  x  10^-^6) ') 
xlabel ( ' \itPv\rm  (GPa  m/s)') 
ylabel ( ' \itk_A\rm' ) 

legend (' \itk_A\rm' ,' \itk_A\rm  adjusted' ,' \itk_A\rm  averaged') 
set (gcf , 'Color', 'w') 

%  Now  recalculate  the  wear  with  the  newly  established  kA_Average 
Wnew=kA  avg ' . *p* 12 *An . /H; 

Wnew  SI=Wnew*53731 . 6; 

figure (3)  %  Thesis  Figure  18 

%plot (pv, Wexp, ' b . ' , pv, Wnew, ' g . ' ) 

plot (pv  SI , Wexp  SI , ' b . ' , pv  SI , Wnew  SI , ' g . ' ) 

%xlabel ( ' \itPv\rm  (psi  ft/s  x  10^-^6) ') 
xlabel (' \itPv\rm  (GPa  m/s)') 

%ylabel  (  '  Wear  Rate  (in''3/ft)') 
ylabel ('Wear  Rate  (mm^3/m) ' ) 
set (gcf , 'Color', 'w') 

legend (' Experimental ' , 'New  \itk_A\rm') 


5-9-5-Q-5-9-9-Q-9-9-5-Q-5-9-5-Q-9-9-5-Q-5-9-9-Q-5-9-5-9-5-9-5-9-9-9-5-9-5-9-9-9-5-9-5-9-5-9-9-Q-9-9-5-9-5-9-5-9-5-9-5-Q-5-9-5-9-9- 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
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%  Looking  at  the  Coefficient  of  Friction 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

xdata=10 : 8000; 

mu  funcl  =  l .  57  04  .  *  (xdata)  . .  22  99; 

C=min  (mu  fund) -0.02; 
mu  func2=mu  funcl-C; 

figure  (4)  %  Thesis  Figure  19 

set (gcf ,  'Color',  'w') 

plot ( Pvtil , mu, ' b . ' , xdata, mu  fund , ' r ' , xdata, mu  func2 , ' g ' ) 
ylim([0  1.4]) 

set (get (gca, ' XLabel ' ) , ' Interpreter ' , ' Latex ' ) 
xlabel ( '$\tilde{P}\tilde{v}$' ) 

ylabel (' Coefficient  of  Friction  (  \it\mu\rm) ') 

legend (' Montgomery ','$ \mu  =  f ( \tilde { P } \tilde {v} )$',.. . 

'$\mu  =  f  ( \tilde  {  P } \tilde  { v} ) '^  *  $  '  ) 
set ( legend, ' Interpreter ' , ' Latex ' ) 

text ( ' Interpreter '  ,  ' Latex '  ,  ' string '  ,  . . . 

'$\mu  =  1 . 57  04  ( \tilde  {  P } \tilde  { v} ) '^  { -0 . 22  9  9  }  $  '  ,  '  position 
[4000, .5, ] ) 
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%  Script  Montgomery_Melt .m 

%  Capt  Greg  Cameron 
%  Thesis  -  Winter  2007 

%  This  script  reads  in  the  Montgomery  data  from  table  10  and 
%  reproduces  the  Lim  and  Ashby  melt  model  and  plots  the  Montgomery 
%  data  on  that  model  and  then  rotates  in  3D  (Thesis  Figure  28) 

g, 

o 

9-9-5-Q-5-9-9-Q-5-9-5-9-5-9-9-Q-5-9-5-Q-5-9-5-Q-5-9-5-9-5-9-5-Q-5-9-Q-9-9-9-5-Q-5-9-5-Q-5-9-9-Q-5-9-5-Q-5-9-5-Q-5-9-5-9-9-9-9-Q-5-9-5-9-5- 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  lb  =  really  affects  how  this  data  matches  Lim  and  Ashby 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


clc 

clear  all 
close  all 


warning  off  all 
format  long 


%Definition  of 
a  =  .014105; 

H  =  1.45e5; 

TO  =  527.4; 
Tstar  =  399.6; 
Tm  =  3240; 

L  =  3.05e5; 


Constants  common  for  both  Lim  and  Ashby 
%Thermal  Diffusivity  (Steel) 

%Steel  Hardness 

%Sink  Temperature  (20  C  Lim/Ashby) 
%222  K  Lim  and  Ashby 
%Melting  Point  for  Steel 
%Latent  heat  of  melting  for  steel 
%  (Lim  and  Ashby  2.1e9  J/m''3)  -  psi 


and  Montgomery 
in''2  /s 
psi 
R 
R 
R 


%f ile= [ ' z : \Greg\Work\Air  Force\AFIT\Thesis\MATLAB\Mont5 . txt ' ] ; 
f ile= ' Month . txt ' ; 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  This  is  for  the  Montgomery  Data 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%velocity  ( f t/s )  ,  pressure  (psi )  ,  ( - )  ,  distance  (ft)  ,wear(in''3/ft*le-7) 

[v, p, mu, d, Wexp]  =  textread (file, ' %f  %f  %f  %f  %f'); 

%Variables  used  by  Montgomery 
r0=0.040; 

An=r0^2*pi; 

lb=0.04; 

beta=6; 

Wexp=Wexp. /An/12*le-7; 

Ftil=p/H; 

vtil=v .*r0./a.*12; 

Q-9-5-Q-5-9-5-Q-9-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-9-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-9-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5- 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Calculate  Melt  wear  from  Montgomery  using  Lim  and  Asbhy  equations 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

alpha=l . / (2  +  lb* (pi/8/a/rO . *v*12)  . ^  ( .  5) ) ; 

Wmelt= ( (Tm-TO) /Tstar*H/L/beta) . /vtil . * (alpha. *mu. *Ftil . *vtil* . . . 
(Tstar*beta/ (Tm-TO) ) -1) *An*12; 

5-9-5-9-5-9-5-Q-5-9-5-Q-9-9-5-Q-5-9-5-Q-9-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-9-Q-5-9-5-Q-5-9-5-Q-5-9-5-Q-9-9-5-Q-5-9-5- 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Generating  Data  to  plot  contour  lines  for  Lim  and  Ashby  Model 

Q-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-9-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5-9-5- 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

rO  model=.04;  %Stated  in  Table  3  Lim  and  Ashby  r0=1.5  mm 

An  model=r0  model^2*pi; 
lb  model= . 004 ; 
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beta_model=l 0 ; 

X=logspace (-1, 5, 100) ;  %This  is  Vtilde 
Y=logspace (-4, 1, 100) ;  %This  is  Ftilde 

[XI, YI]  =  meshgrid (X, Y) ; 

vXY=XI *a* 12 /rO  model; 

muXY=0 . 78-0 . 13*logl0 (XI) ;  %Using  Lim  and  Ashby's  equation  9 
alphaXY=l .  /  (2  +  lb  model*  (pi/8/a/rO  model .  *vXY*  12 )  . '^  ( .  5 )  )  ; 

WmeltXY= ( (Tm-TO) /Tstar*H/L/beta  model) . /XI . * (alphaXY. *muXY. *YI . * . . . 
XI* (Tstar*beta  model/ (Tm-TO) ) -1) *An  model*12; 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Generate  Plot 

9-9-5-9-5-9-9-Q-5-9-5-Q-5-9-5-Q-9-9-5-Q-5-9-5-Q-5-9-5-Q-5-9-9-Q-9-9-5-9-5-9-9-9-5-9-5-Q-5-9-9-Q-9-9-5-Q-5-9-9-Q-5-9-5-Q-5-9-9-Q-9-9-5-9- 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


s=[le-7  le-6  le-5  le-4  le-3] ; 
contours (XI , YI , WmeltXY,  s ) 

set (get (gca,  ' XLabel ' ) ,  ' Interpreter ' ,  ' Latex ' ) 
set (get (gca,  ' XLabel '),  'String',  '$\tilde{v}$') 
set (get (gca, ' YLabel ' ) , ' Interpreter ' , ' Latex ' ) 
set (get (gca, ' YLabel '), 'String', '$\tilde{F}$') 
set  (get (gca,  ' ZLabel ' ) ,  ' Interpreter ' ,  ' Latex ' ) 

set (get (gca,  ' ZLabel '),  'String',  '$\tilde{w}$',  ' Font Size ',16,  ' Rotation ', 0 

) 

zlim([le-7  le-2]) 

set (gca, 'XScale', 'log') 

set (gca,  'XTick',  [le-1  leO  lei  le2  le3  le4  le5]  ) 
set (gca,  'YTick',  [le-4  le-3  le-2  le-1  leO  lei]) 
set (gca, 'YScale', 'log') 
set (gca,  'ZScale',  'log') 

set (gca,  ' XMinorGrid '  ,  ' of f '  ,  ' YMinorGrid '  ,  ' of f '  ,  ' ZMinorGrid ' ,  ' of f ' ) 
set (gca,  ' XMinorTick '  ,  ' of f '  ,  ' YMinorTick  '  ,  ' of f '  ,  ' ZMinorTick ' ,  ' of f  '  ) 
set (gcf , 'Color', 'w') 
hold  on 

%  Now  plot  experimental  Montgomery  Data  as  points 
plots (vt 11, Ft 11, Wexp, ' b . ' ) 
hold  on 

1  =  0; 

J=90; 

set (gca, 'View' , [I  J] ) 

Az=l; 

El=91; 

KAz=.3;  %Gain 

KE1=1;  %Gain 

mov  =  avif lie ( ' Montgomery_Melt . avi  '  ,  ' FPS ' , 2 0 ) ; 
for  1=1:155 

Az=Az-KAz ; 

E1=E1-KE1; 
if  El<50 
KAz=l; 
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KE1=0; 


0, 

o 


end 

mov 


end 

if  Az<-112 

KAz=0; 

KE1=-1; 

end 

Look;=  [Az  El ]  ; 

set (gca,  ' View Look)  ; 

M(i)  =  getframe (gcf ) ; 
mov  =  addf rame (mov, M ( i ) ) ; 

=  close (mov) ; 
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Appendix  B, 


CTH  Input  File 


*eor*  genin 
* 

*  cthgen  input  for  wear  simulation 

* 

* 

*  I  - >  I 

*  I  I  I 

*  I  V  / 

■k  _ 

* 

* 

*  vx=varies,  vy=-l  m/ s  V300  Steel  Slider,  1080  Steel  Rail,  No  Atm. 

*  No  Slide  line.  mix=l  frac=l  Rounded  corner. 

*  Added  mass  on  top  to  simulate  sled  mass 

* 

-1.98  m/s,  798  m/sec,  coating,  asperity,  .002  cm  mesh,  T0=Profile 

control 

mmp3  *  enable  multiple  material  temps  and  pressures  in  each  cell 

ep  *  enable  elastic-plastic  computations 

vpsave  *  save  cell  yield  stress  and  plastic  strain  rate  data 

endcontrol 

*  MESH  AND  BLOCK  DEFINITION  SET 

■k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

*  geom=2DR ( rectangular  x,y) 

*  geom=2DC (cylindrical  x=radius,  y=axis) 

*  geom=3DR ( rectangular  x,y,z) 

*  type=e  (Eulerian) 

*  x#=coordinate  range  for  plot 

*  y#=coordinate  range  for  plot 

*  dxf=width  of  first  cell  in  the  region 

*  dxl=width  of  last  cell  in  the  region 

*  n=number  of  cells  added  in  this  region 

*  w=total  width  of  this  region  in  centimeters 

*  r=ratio  of  adjacent  cell  widths 

"'JJJJJJJJJJJJJJJJJJJJJJJ 

mesh 

block  1  geom=2dr  type=e  *  coord  for  2D  rect  Eulerian  mesh 


xO  = 

xl 

0.0000 
w  =  15.00 

dxf  = 

1 . 000 

dxl 

=  0.1 

x2 

w  =  3.000 

dxf  = 

0.100 

dxl 

=  0.050 

x3 

w  =  2.400 

dxf  = 

0.050 

dxl 

=  0.002 

x4 

w  =  1.500 

dxf  = 

0.002 

dxl 

=  0.002 

x5 

w  =  3.000 

dxf  = 

0.002 

dxl 

=  0.100 

x6 

w  =  5.100 

dxf  = 

0.100 

dxl 

=  1.000 

endx 
yO  = 

-4 . 000 
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yi 

w  =  2.200 

dyf  = 

0.100 

dyl  = 

0.050 

y3 

w  =  1.000 

dyf  = 

0.050 

dyl  = 

0.002 

y4 

w  =  1.600 

dyf  = 

0.002 

dyl  = 

0.002 

y5 

w  =  2.200 

dyf  = 

0.002 

dyl  = 

0.050 

y6 

w  =  5.000 

dyf  = 

0.050 

dyl  = 

0.100 

y7 

w  =  41.50 

dyf  = 

0.100 

dyl  = 

1.000 

endy 

endblock 
endmesh 

insertion  of  material 
block  1 

package  rail 

material  2 
numsub  100 
temperature 
xvel  0 . 0 
yvel  0.0 
insert  box 

pi  0.0  0.0 
p2  50.0  -4.0 
endinsert 
endpackage 

package  coating  *  normal  level  rail  coating 
material  1 
numsub  100 

temperature  =  2.55935e-2  *  eV  =  74.93F  =  297  K 
xvel  0 . 0 
yvel  0.0 
insert  box 
pi  0.0  0.0 
p2  50.0  0.01 
endinsert 
delete  triangle 
pi  20.46  0.0 
p2  24.0  0.0 
p3  20.63  0.03 
enddelete 
endpackage 

package  coating2 
material  1 
numsub  100 

temperature  =  2.55935e-2  *  eV  =  74.93F  =  297  K 
xvel  0.0 
yvel  0.0 
insert  triangle 
pi  20.46  0.01 
p2  24.0  0.01 
p3  20.63  0.04 
endinsert 
delete  triangle 
pi  20.46  0.0 


*  normal  level  rail  height 

=  2.55935e-2  *  eV  =  74.93F  =  297  K 
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p2  24.0  0.0 
p3  20.63  0.03 
enddelete 
endpackage 

*This  part  is  for  the  asperity  part  of  the  rail 
package  rail2 

material  2 
numsub  100 

temperature  =  2.55935e-2  *  eV  =  74.93F  =  297  K 
xvel  0 . 0 
yvel  0.0 

insert  triangle  *This  should  match  the  delete  triangle  above, 
pi  20.46  0.0 
p2  24.0  0.0 
p3  20.63  0.03 
endinsert 
endpackage 

package  shoel  *  shoe  on  normal  level  rail 

material  3 
numsub  100 

temperature  =  6.489e-2  *  eV  =  753.1  K 

xvel  7 . 98e+4 
yvel  -1 . 98e+2 
insert  box 
pi  0.300,  0.01 
p2  20.62,  0.03 
endinsert 
delete  circle 

center  20.42,  0.21 
radius  0.2 
enddelete 
delete  box 

pi  20.42,  0.21 

p2  20.62,  .01 

enddelete 

endpackage 

package  shoe2  *  shoe  on  normal  level  rail 

material  3 
numsub  100 

temperature  =  5.795e-2  *  eV  =  672.4  K 
xvel  7 . 98e+4 
yvel  -1 . 98e+2 
insert  box 
pi  0.300,  0.03 
p2  20.62,  0.045 
endinsert 
delete  circle 

center  20.42,  0.21 
radius  0.2 
enddelete 
delete  box 

pi  20.42,  0.21 


no 


.01 


p2  20.62, 

enddelete 
endpackage 

package  shoe3  *  shoe  on  normal  level  rail 

material  3 
numsub  100 

temperature  =  5.495e-2  *  eV  =  637.7  K 
xvel  7 . 98e+4 
yvel  -1 . 98e+2 
insert  box 
pi  0.300,  0.045 
p2  20.62,  0.075 
endinsert 
delete  circle 

center  20.42,  0.21 
radius  0.2 
enddelete 
delete  box 

pi  20.42,  0.21 

p2  20.62,  .01 

enddelete 

endpackage 

package  shoe4  *  shoe  on  normal  level  rail 

material  3 
numsub  100 

temperature  =  4.927e-2  *  eV  =  571.8.0  K 
xvel  7 . 98e+4 
yvel  -1 . 98e+2 
insert  box 
pi  0.300,  0.075 
p2  20.62,  0.09 
endinsert 
delete  circle 

center  20.42,  0.21 
radius  0.2 
enddelete 
delete  box 

pi  20.42,  0.21 

p2  20.62,  .01 

enddelete 

endpackage 

package  shoe5  *  shoe  on  normal  level  rail 

material  3 
numsub  100 

temperature  =  4.672e-2  *  eV  =  542.2  K 
xvel  7 . 98e+4 
yvel  -1 . 98e+2 
insert  box 
pi  0.300,  0.090 
p2  20.62,  0.135 
endinsert 
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delete  circle 

center  20.42,  0.21 
radius  0.2 
enddelete 
delete  box 

pi  20.42,  0.21 

p2  20.62,  .01 

enddelete 

endpackage 

package  shoe6  *  shoe  on  normal  level  rail 

material  3 
numsub  100 

temperature  =  4.038e-2  *  eV  =  468.6  K 
xvel  7 . 98e+4 
yvel  -1 . 98e+2 
insert  box 
pi  0.300,  0.135 
p2  20.62,  0.165 
endinsert 
delete  circle 

center  20.42,  0.21 
radius  0.2 
enddelete 
delete  box 

pi  20.42,  0.21 

p2  20.62,  .01 

enddelete 

endpackage 

package  shoe7  *  shoe  on  normal  level  rail 

material  3 
numsub  100 

temperature  =  3.712e-2  *  eV  =  430.8  K 
xvel  7 . 98e+4 
yvel  -1 . 98e+2 
insert  box 
pi  0.300,  0.165 
p2  20.62,  0.225 
endinsert 
delete  circle 

center  20.42,  0.21 
radius  0.2 
enddelete 
delete  box 

pi  20.42,  0.21 

p2  20.62,  .01 

enddelete 

endpackage 

package  shoe8  *  shoe  on  normal  level  rail 

material  3 
numsub  100 

temperature  =  3.240e-2  *  eV  =  375.9  K 
xvel  7 . 98e+4 
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yvel  -1 . 98e+2 
insert  box 
pi  0.300,  0.225 
p2  20.62,  0.315 
endinsert 
delete  circle 

center  20.42,  0.21 
radius  0.2 
enddelete 
delete  box 

pi  20.42,  0.21 

p2  20.62,  .01 

enddelete 

endpackage 

package  shoe9  *  shoe  on  normal  level  rail 

material  3 
numsub  100 

temperature  =  2.842e-2  *  eV  =  329.8  K 
xvel  7 . 98e+4 
yvel  -1 . 98e+2 
insert  box 
pi  0.300,  0.315 
p2  20.62,  0.750 
endinsert 
delete  circle 

center  20.42,  0.21 
radius  0.2 
enddelete 
delete  box 

pi  20.42,  0.21 

p2  20.62,  .01 

enddelete 

endpackage 

package  shoelO  *  shoe  on  normal  level  rail 

material  3 
numsub  100 

temperature  =  2.55935e-2  *  eV  =  297  K 
xvel  7 . 98e+4 
yvel  -1 . 98e+2 
insert  box 
pi  0.300,  0.750 
p2  20.62,  2.550 
endinsert 
delete  circle 

center  20.42,  0.21 
radius  0.2 
enddelete 
delete  box 

pi  20.42,  0.21 

p2  20.62,  .01 

enddelete 

endpackage 
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package  shoetip 
material  3 
numsub  100 

temperature  =  2.55935e-2  *  eV  =  74.93F  =  297  K 
xvel  7 . 98e+4 
yvel  -1 . 98e+2 
insert  circle 

center  20.42  0.21 
radius  0.2 
endinsert 
endpackage 

package  sledsim  *  sled  mass  concentrated  on  top  of  shoe 

material  4 
numsub  100 

temperature  =  2.55935e-2  *  eV  =  74.93F  =  297  K 
xvel  7 . 98e+4 
yvel  -1 . 98e+2 
insert  box 

pi  0.300,  49.08 
p2  20.62,  2.55 
endinsert 
endpackage 

endblock 

endinsertion 

edit 

block  1 
expanded 
endblock 
endedit 

*  TRACER  DEFINITION  SET 

tracer  4x4  set  of  points  in  the  bottom  of  the  shoe 
*Bottom  Row  in  shoe 

add  2.84,  .010  to  18.08,  .010  n=7 
add  2.84,  .020  to  18.08,  .020  n=7 
*2nd  Row 

add  2.84,  .060  to  18.08,  .060  n=7 

*3rd  Row 

add  2.84,  .180  to  18.08,  .180  n=7 

*4th  Row 

Add  2.84,  .540  to  18.08,  .540  n=7 

*5th  Row 

Add  2.84,  40.00  to  18.08,  40.00  n=7 

endt 

*  EQUATION  OF  STATE  DEFINITION  SET 

"'JJJJJJJJJJJJJJJJJ 

eos  *  Sesame  tabular  option  (primary  MIG  model) 
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*  MATl  SES  IRON 

MATl  SES  GREPXYl  *for  epoxy 
MAT2  SES  IRON 
MATS  SES  STEEL_V300 
MAT4  MGR  PLATINUM 
endeos 


epdata 
vpsave 
Istrain 
mix  =  3 


*  elastic-plastic  data 

*  cell  yield  stress  and  plastic  strain  rate  data  is  saved 

*  compute  and  save  Lagrangian  strain  tensor  components 

*  vol  avg  yield  strength  normalized  by  sum  of  vol  fractions 


matep  = 

1  ■> 

'Epoxy 

Glider  Coating 

poisson  0.46 

yield 

l.OeS 

matep  = 

2  ZE  USER 

■k 

1080  Steel  rail 

aze  = 

0.825el0 

k 

dynes/cm' 

'2 

nze  = 

0.289 

cl  ze  = 

^  4.0el0 

k 

dynes/cm' 

'2 

c2  ze  = 

^  0.0 

c3ze  = 

^  160.0  * 

1/eV 

c4  ze  = 

^  12.0 

k 

1/eV 

c5ze  = 

^  0.266el0 

k 

dynes/cm' 

'2 

poisson  0.27 

matep  = 

3  ZE  USER 

k 

VascoMax 

300  shoe 

aze  = 

1.42el0 

k 

dynes/cm' 

'2 

nze  = 

0.289 

cl  ze  = 

^  4.0el0 

k 

dynes/cm' 

'2 

c2  ze  = 

^  0.0 

c3ze  = 

^  79.0 

k 

1/eV 

c4  ze  = 

^  3.0 

k 

1  /ev 

c5ze  = 

^  0.266el0 

k 

dynes/cm' 

'2 

poisson  0.283 

matep  = 

4  -k 

what  . 

is  this  for? 

poisson  .2 

yield 

lOelO 

* 

* 

*  BL  algorithm  inputs 


*  blint  1 

*  soft  1 

*  hard  3 

*  csl  1.5 

*  cbl  1.5 

*  f ric  0 . 3 

*  corr 

*  nofreeze 

*  ende 

ende  *  end  elastic-plastic  data 
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*eor*  cthin 

-1.0  m/s,  100.0  m/sec,  coating,  .002  cm  mesh 

"'JJJJJJJJJJJJJJJJJJJJJJJJJJJJ 

*C0NTR0L  DATA  SET 

'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

*  mmp  enables  multiple  material  temps 

*  and  pressures  in  mixed  cells 

*  frac=l  changes  fracture  default 

'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

*  CELL  THERMODYNAMICS  INPUT  SET 

'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

*  cellthermo 

*  dtmax  =  max  temp  difference  allowed  in  mixed  cells 

*  mmp  distributes  volume  and  energy  based  on  volume 

*  fractions  of  material  in  the  cell 

*  mmp=default 

*  mmpl=same  as  mmp  except  uses  new  logic 

*  and  distributes  volume  and  energy 

*  based  on  volume  fraction  cubed  divided 

*  by  mass  of  material  in  the  cell 

*  mmp2=allocation  of  work  done  on  the  cell  is 

*  dependent  on  material  compressibility 

*  and  allows  pressure  relaxation  between  materials 

*  in  a  cell 

*  endc 

*//////////////////////////// 


control 

mmp3  *  allocates  PdV  work  done  between  materials  of  a  mixed  cell 

tstop  =  lO.OOE-6  *  stop  criteria  -  total  simulation  time 

nscycle  =  55000  *  maximum  number  of  cycles  to  be  run 
rdumpf  =  991800.  *  time  for  back-ups  of  restart  file  updates 

cpshift  =  600.  *  time  left  when  normal  termination  procedure  begins 

ntbad  =  99999999 

dtcourant  =0.6  *  Courant  condition  multiplier 
ygravity=-980  *  Acceleration  due  to  gravity 

endc 

Convct 

convection 
residual 
interface 
endc 

fracts  *  enable  fracture  data  (dynes/cm^2) 

*pfracl  =  -0.5e7  *  fracture  stress  or  pressure  for  nth  material 

pfracl  =  -1.0e8 

pfrac2  =  -2.0el0 

pfrac3  =  -2.5el0 

pfrac4  =  -1.2el0 

pfmix  =  -1.20el0  *  fract  stress  or  pressure  in  a  cell  with  no  voids 


*  enable  convection  of  internal  energy 
=1  *  use  int  energy  slope  and  mass  density,  discard  KE 

=  high  resolution  *  scheme  for  interface  tracker 
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pfvoid  =  -1.20el0  *  fract  stress  or  pressure  in  a  cell  with  a  voids 
endf racts 

edit 

shortta  *  short  edits  based  on  time 

time  =0.0  ,  dt  =  l.Oe-6 

ends 

longt  *  long  edits  based  on  time 

time  =  O.OeO  ,  dt  =  1.0 
endl 

plott  *  plot  dumps  based  on  time 

time  O.Oe-6  dtfrequency  0.5e-6 
endp 

histt  *  tracer  history  based  on  time 

time  O.Oe-6  dtfrequency  0.5e-6 
htracer  all 
endhistt 
ende 

"'JJJJJJJJJJJJJJJJJJJJJJJJJJJJ 

*  Hydrodynamic  Boundary  Conditions 

'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

*  0=symmetry 

*  l=sound  speed  based  absorbing 

*  2=extrapolated  pressure  with  no  mass  allowed  to  enter 

*  3=extrapolated  pressure  but  mass  is  allowed  to  enter 

J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J  J 

boundary  *  enable  boundary  condition  data 

bhydro  *  enable  hydrodynamic  boundary  conditions 

block  1 

bxbot  =  1  ,  bxtop  =1  *  X  =  xmin,  x  =  xmax,  1  =  sound  speed 

based  absorbing  (transmitting)  b.c. 

bybot  =  1  ,  bytop  =1  *  y  =  ymin,  y  =  ymax,  1  =  sound  speed 

based  absorbing  (transmitting)  b.c. 
endb 
endh 
endb 

*heat conduct ion 

*  MATl  TABLE  =  3 

*  MAT2  TABLE  =  1 

*  MAT 3  TABLE  =  2 

*endh 

*  DEFTABLE=1  *  1080  STEEL 

*T (eV)  k (erg/s/eV/cm) 

*  1.4684e-3  4.7700el0 

*  1.0377e-2  4.8100el0 

*  1.9090e-2  4.5200el0 

*  2.7900e-2  4.1300el0 

*  3.6711e-2  3.8100el0 

*  4.5521e-2  3.5100el0 

*  5.4332e-2  3.2700el0 

*  6.3142e-2  3.0100el0 


*  enable  heat  conduction 

*  condty  tabls  defined  in  DEFTABLE  below 
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*  7.1953e-2  2.4400el0 

*  8.9574e-2  2.6800el0 

*  l.lllle-1  S.OlOOelO 

*  endd 

*  DEFTABLE=2  *  VascoMax  300  Steel 

*T (eV)  k (erg/s/eV/cm) 

*  3.6711e-3  2.4715el0 

*  1.4684e-2  2.7424el0 

*  2.9369e-2  2.9794el0 

*  3.9158e-2  3.0132el0 

*  endd 

*  DEFTABLE=3  *  Epoxy 

*T (eV)  k (erg/s/eV/cm) 

*  3.6711e-3  6.5e8 

*  1.4684e-2  6.5e8 

*  2.9369e-2  6.5e8 

*  3.9158e-2  6.5e8 

*  endd 


*vadd  *  Added  velocity  to  maintain  gouging  in  view 

*  block  =  1 

*  tadd  =0.0 

*  xvel  =  -1.08333e+5 
*endvadd 

*mindt  *  minimum  allowable  time  step  in  mesh 

*  time  =  0.0  dt  =  l.e-12 
*endm 

maxdt  *  maximum  allowable  time  step  in  mesh 

time  =  0.0  dt  =  .01 
endm 
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Appendix  C, 


Calculation  of  Sled  Mass  for  CTH 


Evenly  distribute  the  mass  over  four  shoes  and  divide  by  the  width  to  make  it  plane  strain. 


Length  of  the  Shoe: 


1854 

mass  := - lb 

4-(10.1Q 


mass  =  45.621b 


Length  :=  S  in 


Height  and  Density  of  the  Respeetive  Shoe  (VaseoMax)  and  the  Mass  (platinum): 


HeighUm„x:=2-54cm 


Heightpiatinum  :=  (49.07-  2.54)-c 


P  platinum  ■ 


:=21.45- 


Caleulation  of  the  volume  for  eaeh  material: 


Volvmax:=  LengthHeight^j^^^-  1cm 


Volpiatinum  :=  Length  Heightpi^tinum'  Lt 


VoVmax=5L61cm 


Volpiatinum  =  945.49cm 


Calculation  of  the  mass  for  each  material: 


Mass  ymax  ■  P  vmax^°Vn 


Mass  platinum  ■  P  platinum' ^°*platir 


Mass  y^ax=  0-91  lb 


Mass  platinum  =  44-71  lb 


Total  mass  should  equal  the  mass  above: 


Mass  :=  Mass  platinum 


+  Mass , 


Mass  =  45.621b 


Therefore  the  platinum  mass  in  CTH  should  have  a  height  of: 


HeigLtpiatinum  =  46.53cm 
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Appendix  D, 


DADS  Data  File 


t 

(s) 

V 

(in/s) 

(Ibf) 

Vv 

(in/s) 

0.0000 

0.000 

0.0 

0.000 

0.0001 

0.000 

0.0 

-0.625 

0.0002 

0.000 

0.0 

-2.945 

0.0003 

0.006 

0.0 

-5.193 

0.0004 

0.125 

0.0 

-7.082 

0.0005 

1.079 

0.0 

-8.784 

0.0006 

5.450 

0.0 

-9.279 

0.0007 

18.651 

0.0 

-9.372 

0.0008 

46.947 

0.0 

-9.752 

0.0009 

91.282 

0.0 

-10.558 

0.0010 

141.067 

0.0 

-12.785 

0.0011 

175.615 

0.0 

-15.607 

0.0012 

176.052 

0.0 

-17.325 

0.0013 

140.216 

0.0 

-18.504 

0.0014 

87.804 

0.0 

-19.287 

0.0015 

49.531 

0.0 

-19.321 

0.0016 

46.915 

0.0 

-19.352 

0.0017 

77.466 

0.0 

-19.971 

0.0018 

116.671 

0.0 

-20.866 

0.0019 

135.712 

0.0 

-21.695 

0.0020 

121.851 

0.0 

-21.972 

0.0021 

86.434 

0.0 

-21.934 

0.0022 

55.443 

0.0 

-21.625 

0.0023 

50.871 

0.0 

-20.865 

0.0024 

76.668 

0.0 

-19.530 

0.0025 

117.567 

0.0 

-17.619 

0.0026 

149.551 

0.0 

-15.311 

0.0027 

154.558 

0.0 

-12.933 

0.0028 

131.215 

0.0 

-10.826 

0.0029 

95.965 

0.0 

-9.204 

0.0030 

73.676 

0.0 

-8.111 

0.0031 

82.509 

0.0 

-7.442 

0.0032 

121.923 

0.0 

-6.947 

0.0033 

171.855 

0.0 

-6.275 

0.0034 

204.855 

0.0 

-5.092 

0.0035 

204.214 

0.0 

-3.232 

0.0036 

175.786 

0.0 

-0.777 

0.0037 

144.323 

0.0 

1.954 

0.0038 

136.187 

0.0 

4.473 

0.0039 

161.043 

0.0 

6.273 

0.0040 

206.057 

0.0 

7.036 

0.0041 

246.029 

0.0 

4.628 

0.0042 

260.890 

0.0 

2.468 

0.0043 

247.897 

0.0 

3.191 

0.0044 

221.513 

0.0 

4.820 

0.0045 

203.067 

0.0 

6.973 

0.0046 

207.946 

0.0 

9.042 

0.0047 

237.799 

0.0 

10.431 

0.0048 

280.773 

0.0 

10.779 

0.0049 

318.210 

0.0 

5.841 

0.0050 

333.969 

0.0 

5.529 

0.0051 

322.464 

0.0 

6.584 

0.0052 

291.953 

0.0 

8.062 

0.0053 

260.661 

0.0 

9.328 

0.0054 

246.651 

0.0 

9.828 

0.0055 

257.018 

0.0 

5.575 

0.0056 

283.772 

0.0 

3.344 

0.0057 

309.638 

0.0 

4.333 

0.0058 

319.710 

0.0 

6.291 

0.0059 

310.713 

0.0 

8.685 

0.0060 

291.846 

0.0 

10.806 

0.0061 

277.594 

0.0 

12.075 

0.0062 

278.225 

0.0 

11.717 

0.0063 

294.278 

0.0 

6.525 

0.0064 

317.589 

0.0 

6.859 

0.0065 

336.832 

0.0 

8.146 

0.0066 

343.657 

0.0 

9.873 

0.0067 

336.577 

0.0 

11.354 

0.0068 

321.506 

0.0 

11.963 

0.0069 

308.690 

0.0 

6.987 

0.0070 

306.927 

0.0 

4.071 

0.0071 

317.990 

0.0 

4.993 

0.0072 

335.282 

0.0 

6.816 

0.0073 

348.463 

0.0 

8.967 

0.0074 

351.145 

0.0 

10.725 

0.0075 

345.849 

0.0 

11.437 

0.0076 

342.238 

0.0 

4.160 

0.0077 

349.871 

0.0 

1.037 

0.0078 

371.078 

0.0 

1.966 

0.0079 

399.384 

0.0 

3.436 

0.0080 

424.373 

0.0 

4.720 

0.0081 

438.761 

0.0 

5.179 

0.0082 

442.229 

0.0 

-1.609 

0.0083 

439.991 

0.0 

-3.721 

0.0084 

438.439 

0.0 

-2.983 

0.0085 

441.511 

0.0 

-2.057 

0.0086 

449.479 

0.0 

-1.645 

0.0087 

459.631 

0.0 

-7.772 

0.0088 

467.885 

0.0 

-9.850 

0.0089 

470.843 

0.0 

-8.994 

0.0090 

467.819 

0.0 

-7.752 

0.0091 

461.698 

0.0 

-6.791 

0.0092 

457.655 

0.0 

-7.898 

0.0093 

460.160 

0.0 

-16.952 

0.0094 

470.282 

0.0 

-17.284 

0.0095 

485.381 

0.0 

-16.158 

0.0096 

501.237 

0.0 

-14.731 

0.0097 

514.466 

0.0 

-13.626 

0.0098 

523.070 

0.0 

-13.544 

0.0099 

525.474 

0.0 

-19.808 

0.0100 

520.398 

0.0 

-20.530 

0.0101 

508.769 

0.0 

-19.718 

0.0102 

495.621 

0.0 

-18.721 

0.0103 

488.730 

0.0 

-18.089 

0.0104 

493.523 

0.0 

-20.347 

0.0105 

508.083 

0.0 

-25.198 

0.0106 

523.176 

0.0 

-24.978 

0.0107 

528.454 

0.0 

-24.445 

0.0108 

520.529 

0.0 

-24.104 

0.0109 

506.386 

0.0 

-26.083 

0.0110 

498.940 

0.0 

-28.857 

0.0111 

507.886 

13.6 

-28.315 

0.0112 

532.735 

1199.9 

-24.776 

0.0113 

563.032 

2813.6 

-13.821 

0.0114 

585.157 

4050.4 

-16.217 

0.0115 

590.708 

6087.3 

-14.299 

0.0116 

581.277 

4032.5 

2.319 

0.0117 

567.324 

4791.9 

-0.237 

0.0118 

562.175 

3369.9 

5.566 

0.0119 

574.274 

2558.1 

11.606 

0.0120 

601.699 

2270.9 

8.568 

0.0121 

632.432 

1701.2 

10.488 

0.0122 

651.183 

1142.5 

12.334 

0.0123 

649.405 

881.2 

9.213 

0.0124 

631.933 

668.7 

10.730 

0.0125 

614.756 

416.6 

11.194 

0.0126 

614.304 

204.6 

7.984 

0.0127 

635.519 

45.7 

8.184 

0.0128 

667.838 

0.0 

7.669 

0.0129 

692.652 

0.0 

3.501 

0.0130 

696.747 

0.0 

2.257 

0.0131 

681.517 

0.0 

1.872 

0.0132 

661.327 

0.0 

-1.909 

0.0133 

652.983 

0.0 

-3.738 

0.0134 

664.558 

0.0 

-3.852 

0.0135 

691.108 

0.0 

-6.220 
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Appendix  E 


Montgomery  Data  [19] 


V 

P 

d 

w 

(ft/s) 

(psi) 

(-) 

(ft) 

(in^/ft  X  10’) 

50 

11300 

0.26 

5.7 

200 

600 

11900 

0.27 

7.2 

6 

900 

11500 

0.22 

6.1 

13 

1200 

10900 

0.25 

8.3 

23 

1200 

10500 

0.3 

4.6 

52 

150 

11700 

0.35 

1.5 

250 

900 

11500 

0.24 

10.2 

14 

1200 

11100 

0.27 

13.6 

47 

600 

11800 

0.24 

6.1 

18 

600 

11900 

0.23 

16.2 

9 

600 

11900 

0.23 

3.7 

13 

600 

11100 

0.26 

11.2 

15 

150 

10900 

0.42 

3 

250 

900 

10700 

0.35 

13 

71 

1200 

10700 

0.3 

19.6 

72 

900 

11900 

- 

24 

48 

900 

11500 

0.24 

21.9 

22 

900 

11900 

- 

18.4 

19 

1200 

10900 

0.29 

26.6 

75 

900 

11100 

0.24 

21 

19 

30 

11900 

0.48 

0.36 

86 

60 

10100 

0.63 

0.52 

750 

60 

10900 

0.54 

0.64 

920 

10 

12700 

0.49 

0.081 

96 

20 

12500 

0.35 

0.22 

no 

150 

6000 

0.52 

5.6 

150 

30 

3900 

1.3 

5.7 

160 

30 

3900 

1.3 

2.7 

220 

300 

3900 

0.5 

29 

8 

300 

3700 

0.46 

13.2 

8 

450 

16900 

0.32 

9.8 

79 

600 

12900 

0.28 

13.1 

36 

600 

12000 

0.28 

14.2 

27 

600 

5700 

0.24 

7.6 

6 

600 

12100 

0.3 

14.5 

26 

600 

12400 

0.27 

6.4 

18 

600 

3600 

0.28 

27.2 

5 

600 

3500 

0.29 

15.5 

5 

900 

8300 

0.3 

20.7 

25 

1200 

5400 

0.37 

26.5 

27 

1200 

11400 

0.3 

32.2 

83 
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V 

P 

d 

w 

(ft/s) 

(psi) 

(-) 

(ft) 

(in^/ft  X  10’) 

1200 

19400 

0.24 

5.4 

95 

1200 

18200 

0.24 

10.6 

106 

1200 

18100 

0.21 

16.3 

89 

1200 

19700 

0.22 

19.8 

114 

1200 

18700 

0.23 

25.6 

111 

1200 

23500 

0.22 

25.4 

305 

1200 

4400 

0.39 

32.2 

30 

1500 

17500 

0.2 

26.1 

24 

1500 

3400 

0.28 

38.2 

28 

1800 

4000 

0.3 

42.8 

19 

1800 

4800 

0.29 

40 

38 

1800 

17500 

0.19 

18.7 

29 
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Appendix  F. 


Rail  Undulation 


This  appendix  documents  an  effort  that  was  initially  pursued  in  order  to  analyze 
the  dynamics  of  the  sled  as  it  encounters  the  undulation  in  the  rail.  This  information  was 
not  used  due  to  the  fact  that  DADS  provides  the  vertical  velocity.  The  analysis  was 
simplified  by  considering  a  particle  of  constant  mass  traveling  along  a  curved  path. 
Figure  64  below  is  a  depiction  of  a  particle  as  it  travels  along  a  curved  path  from  point  A 
to  point  B.  At  point  yf,  the  particle  has  a  velocity  of  va.  Similarly,  the  particle  has 
velocity  vg  at  point  B. 


From  dynamics,  it  is  know  that  a  particle  traveling  along  a  curved  path  experiences 
accelerations  in  the  form  of 


^  ^  ^  dv 

a  =  a,  +a„  = - 1 - 

dt  r 


V  is  the  tangential  velocity  provided  by  DADS.  Therefore,  the  radius  of  the  rail 


undulation  can  be  solved  for  using  Newton’s  Second  Law  as 


(29) 
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r  = 


(30) 


mv^ 

F 


The  geometric  relations  associated  with  a  particle  traveling  on  a  curve  path  can  be  seen  in 
Figure  64.  This  shows  that  the  angle  between  the  two  points  A  and  B  can  be  solved  for 
by 


e 


Ad 


(31) 


r 

Ad  is  the  distance  the  sled  travels  between  time  steps  and  is  calculated  be  integrating  the 
velocity  time  curve  accomplished  earlier.  Points  A  and  B  can  be  any  two  successive 
points  in  the  DADS  data.  With  the  angle  6,  the  vertical  velocity  can  be  calculated  by 
assuming  that  the  curved  path  is  comprised  of  short  line  segments,  and  that  the  tangential 
velocity  provided  by  DADS  for  the  previous  time  step  will  be  the  velocity  at  which  the 
sled  encounters  the  next  change  in  track  slope.  Resolving  the  components  of  this 
incoming  velocity  vector  gives  the  geometric  relations  for  the  vertical  velocity  as 
depicted  in  Figure  65  and  equation  (32): 


Figure  65  -  Shoe  as  It  Encounters  the  Next  Segment  of  a  Curved  Rail 

(32) 

This  set  of  calculations  can  be  performed  for  the  entire  set  of  DADS  data  in  order 
to  determine  the  vertical  velocity  for  each  time  step.  This  was  carried  out  and  compared 
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with  the  DADS  provided  values  for  vertieal  veloeity.  This  method  resulted  in  vertieal 
veloeities  less  than  what  DADS  ealeulates.  This  suggests  that  this  teehnique  over 
simplifies  the  problem  in  assuming  that  the  mass  is  eonstant  throughout  the  entire  sled 
run  in  addition  to  negleeting  other  system  forees  sueh  as  springs  and  dampers  used  by 
DADS. 
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